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SUMMARY 


Time-dependent  properties  of  uniaxially  fiber  reinforced  mate¬ 
rials  composed  of  linear  viscoelastic  matrix  and  elastic  fibers  are 
investigated.  Sample  calculations  are  given  for  the  static  and  dynamic 
properties  of  a  viscoelastic  fiber  reinforced  material.  In  addition, 
the  behavior  of  fiber  reinforced  viscoelastic  structures  is  investi¬ 
gated,  and  a  number  of  practical  problems  for  beams,  plates,  and  shells, 
subjected  to  static  and  dynamic  loadings,  are  analyzed. 
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1.  INTRODUCTION 


The  purpose  of  the  present  investigation  is  to  study  the  dy¬ 
namic  behavior  of  uniaxially  fiber  reinforced  material  (FRM)  for  the 
case  of  elastic  fibers  and  linear  viscoelastic  matiix. 

A  large  number  of  investigations  have  dealt  with  the  static 
elastic  behavior  of  uniaxially  FRM.  Of  these,  the  most  fruitful  ap¬ 
proach,  in  the  writer's  opinion,  has  been  the  analysis  given  by  Hashin 
and  Rosen  [1],  which  was  based  on  a  geometrical  model  named  the  compos¬ 
ite  cylinder  assemblage  (CCA) .  The  advantages  of  this  approach  are  here 
briefly  summarized: 

(a)  The  geometrical  model  is  based  on  a  random  (though 
special)  geometry. 

(b)  The  mathematical  analysis  of  the  model  is  rigorous. 

(c)  The  analysis  yields  simple  closed  expressions  for 
four  of  the  five  elastic  moduli  of  the  CCA. 

(d)  Experimental  verification  of  the  results  has  been 
excellent. 

The  relative  simplicity  of  the  approach  makes  it  possible  to 
extend  the  approach  to  more  complicated  cases.  Thus,  the  CCA  model  has 
been  used  successfully  by  Hashin  et  al  [2],  to  analyze  the  elastic  prop¬ 
erties  of  biaxially  FRM.  Another  extension  which  is  of  primary  impor¬ 
tance  for  present  purposes  is  a  recent  analysis  by  Hashin  [3]  of  the 
static  viscoelastic  properties  of  a  uniaxially  FRM  on  the  basis  of  the 
CCA  model.  Our  purpose  here  is  to  use  the  same  model  for  prediction  of 
the  dynanuc  properties  of  a  viscoelastic  uniaxially  FRM.  Such  a  study 
is  of  primary  practical  importance  since: 

(a)  In  many  cases,  the  matrix  of  a  FRM  is  time  dependent 
and,  to  a  good  approximation,  linearly  viscoelastic. 

(b)  Frequently,  structural  elements  of  FRM  are  exposed  to 
vibrations . 

Vibration  analysis  on  the  basis  of  elastic  properties  will 
predict  resonant  frequencies  that  should  be  avoided.  However,  the  time- 
dependent  and  dissipative  nature  of  the  matrix  produces  damping  of  vi¬ 
brations,  which  is  of  beneficial  effect.  Free  vibrations  will  rapidly 
die  out,  while  forced  vibrations  will  have  a  less  dangerous  character 
than  that  predicted  on  the  basis  of  pure  elastic  behavior.  Thus,  vibra¬ 
tion  analysis  which  takes  into  account  the  viscoelastic  nature  of  the 
materials  should  lead  to  more  realistic  design  criteria. 


1 


2. 


SUMMARY  OF  STATIC  VISCOELASTIC  BEHAVIOR  OF  UN I AX  I ALL Y  FRM 


The  present  summary  is  based  on  the  results  obtained  in  Refer¬ 
ence  [3].  The  static  stress-strain  relations  of  any  anisotropic  hetero¬ 
geneous  viscoelastic  material  may  be  written  in  the  form 


°u  (t)  ■ 


*  d£kt(T) 

CijkH  (t  -  ’>  -3T-  • 


£ij  (t)  = 


*  ,  ,  dok£(T) 

Sijk£  (t  t)  dT 


(1) 

(2) 


where  o^j  are  average  stresses,  are  average  strains,  C*.^  are 
the  effective  relaxation  moduli,  and  S*j^£  are  the  effective  creep  com¬ 
pliances.  The  range  of  subscripts  is  1,  2,  3,  and  repeated  subscripts 
indicated  summation  over  1,  2,  3.  It  should  be  borne  in  mind  that  be- 
cause  of  the  dual  nature  of  (1)  and  (2),  the  and  Sijkn  are  relat¬ 

ed.  This  relationship  will  be  given  later.  It  is  convenient  to  take 
the  Laplace  transform  (LT)  of  (1)  and  (2).  We  define  the  LT,  <f(p)  of  a 
function  (t)  by 


4>(p)  = 


.00 

e  <Ht)  dt. 
o 


(3) 


By  use  of  the  convolution  theorem  and  the  rules  for  LT  of  derivatives, 
the  LT  of  (1)  and  (2)  assumes  the  form 


5u  (p>  ■  p  <p)  Ek«  (p}  •  (4) 

hi  (p>  ■  p  Vt  <p>  5u  (p)  •  (5) 

It  is  now  seen  that  the  p  6*^0  and  p  S*jk£  matrices  must  be  reciprocal. 
This  is  the  relation  between  tne  effective  relaxation  moduli  and  the 
creep  compliances.  This  relationship  can  be  transformed  to  the  time 
domain.  However,  it  then  becomes  very  complicated  and  is  not  very  use¬ 
ful  in  that  form. 


It  is  convenient  to  define  the  following  quantities: 


* 

'ijk£ 

(p>  '  p  c*jkt  <p)  • 

(6) 

* 

.* 

(7) 

ijk£ 

<p>  ■  p  sijk«  <p)  • 

2 


Then  (4)  and  (5)  assume  the  form 


Vp)  ■  B*jki<p)  Ju<p)  • 

00 

rij(p)  ■  J*jkl(p)  °kt(p)  • 

(9) 

The  relationships  (8)  and  (9)  bear  a  formal  resemblance  to 
elasticity  stress-strain  laws.  Therefore,  the^quantities  B^.,^  are 
called  transform  domain  (TD)  moduli,  and  the  Jijk£,  are  called  TD  compli¬ 
ances. 


The  results  given  so  far  apply  for  any  anisotropic  viscoelas¬ 
tic  material.  A  uniaxially  fiber  reinforced  material  is  transversely 
isotropic,  with  respect  to  an  axis  in  fiber  direction.  As  a  consequence 
of  this  symmetry,  it  may  be  shown  that,  as  in  the  elastic  case,  there 
are  only  five  independent  material  functions  (i.e.,  relaxation  moduli  or 
creep  compliances).  The  relation  (8)  can  then  be  simplified  to  the  form 


ou(p) 

°22<p) 

o33(p) 

o12(p) 


Bll(p)  ^ll(p)  +  B12(p)  ?22(p) 

it  —  it  2 

Bl2(p)  Gn(p)  +  B22^  g22^ 

b12(p)  Gn<P)  +  B23^P)  ^22(p) 

2B44(p)  f12(p)  » 


+  b*2 (p)  c33(p)  , 
*  2 

B23^p)  e33(p)  » 

*  z 

+  B22(p)  c33(p)  , 


(10) 

(ID 

(12) 

(13) 


°13(p) 


2B44<p)  Gn<P) 


13' 


(14) 


o23(p) 


B22^p^  ~  B23 ^  g23^p^  * 


(15) 


where  x^  is  fiber  direction  and  x2,  x3  are  transverse  directions 
(Figure  1).  A  similar  simplification  holds  for  (9).  It  is  seen  that 

it  it 

(10)  through  (15)  contain  the  five  independent  quantities  B..,  B  „, 

*  *  *  11 
B22*  B23’  and  B44‘  In  view  (B)>  we  may  deflne  the  quantities  C.^, 

^12’  ^22’  ^23’  anc*  ^44’  whlch  are  related  to  the  former  quantities  by 


* 

B 


mn 


(P) 


(P) 


mn  =  11,  22,  23,  44.  (16) 
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Figure  1 


X| 


X2 


Geometrical  Model  of  the  Material. 


4 


When  (16)  is  inserted  into  (10)  through  (15)  and  inverted  to  the  time 
domain,  we  have 


5u(t)  = 


'  de  (x)  .  de„_  (x )  .  de_„  (x) 

cn<‘->  — Jt  +  -57—  +  C12<^>  -5T- 


dx  , 
(17) 


S22(t> 


de  (x)  .  de99 (t)  .  de_„(x) 

—57 —  +  c22(t->  “ST—  +  c23‘m>  -57- 


dt  , 
(18) 


S23(t) 


de  (x)  .  de„_ (t )  .  de„  (x) 

C!2<t-)  “ST—  +  C23<'-^  “IT—  +  C22<'->  -E — 


di  , 
(19) 


0 12  <  c “  2 


*  dE12<T) 

-57— 


dt  , 


(20) 


o13(t)  =  2 


*  dE  3(t) 

C44(t->  5t 


dt  , 


(21) 


°23 (t) 


t  r 


o  ^ 


C*2(t-T)  -  C23(t-T) 


d?23(T) 

dx 


dx  . 


(22) 


We  can  write  similar  relations  for  expressions  of  strains  in  terms  of 
stresses.  Thus,  for  example,  the  analogue  of  (17)  would  be 


ft 


Eu(t)  - 


*  do..(x)  .  do.,(x)  .  do--(x) 

“ST—  +  S12<C-’>  “ST—  +  »12l*-T>  47— 


fdT  * 
(23) 


etc. 

The  relations  of  strains  in  terms  of  stresses  involve  the  five  effective 

A  A  A  4r 

creep  compliances  Sll’  S12’  S22’  S23 ’  S44  ,  which  are  related  to  the 
five  relaxation  moduli. 

The  above  choice  of  effective  relaxation  moduli  and  creep  com¬ 
pliances  is  not  very  advantageous  from  a  physical  point  of  view.  To 
bring  out  physically  significant  quantities,  it  is  necessary  to  use 
special  states  of  average  strain  and  average  stress.  The  simplest  sit¬ 
uations  are  presented  by  (20),  (21),  and  (22),  which  all  correspond  to 
pure  shears.  We  adopt  the  notation 


5 


c 

c 


*4(t,  =  G*(t)  ,  (24) 

*2(t)  -  C*3(t)  =  2G*(t)  .  (25) 


It  is  then  seen  that  (20)  and  (21)  are  the  stress  histories  produced  by 
shearing  strain  histories,  where  the  shears  are  in  planes  parallel  to 
the  fibers.  Thus,  Gg(t)  is  the  effective  axial  relaxation  shear 
modulus . 


By  similar  reasoning,  Gx(t)  is  the  effective  transverse  relax¬ 
ation  shear  modulus,  for  shear  in  planes  normal  to  the  fibers. 


The  creep  compliances  ga(t)  and  g£(t)  correspond  to  (24)  and 
(25).  The  meaning  of  these  is  brought  out  by  the  stress-strain  rela¬ 
tions 


^(t)  =  y 


*  d512(T) 

Sa(t->  “if-  dT  • 


(26) 


C13(t)  2 


*  do  (x) 


dx  , 


(27) 


c23(t)  =  2 


do23(T) 


gT(t-T)  . dx  . 


(28) 


Note  that  the  Laplace  transforms  of  (20)  through  (22),  (24),  (25),  and 
(26)  through  (28)  yield  the  relations 


-* 

pGfl(p)  = 


P  8a(p) 


(29) 


pGT(p)  =  — -  .  (30) 

P  8t(P> 

We  now  turn  to  the  more  complicated  relations,  (17)  through 
(19),  between  normal  average  stresses  and  strains.  First,  assume  the 
average  state  of  strain  to  be 

c22(t)  =  t33(t)  =  c(t)  , 

cn(t)  =  c12(t)  =  en(t)  =  e23(t)  =  0  .  (31) 
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This  is  an  isotropic  plane  strain  system,  in  planes  normal  to  the  fibers. 
Introducing  (31)  into  (17)  through  (19),  we  have 


o(t) 

2K^(t) 


o33(t)  =  o ( t )  , 


=  2 


K*(t-x)  d^T-)  di  , 
T  dx  * 


c*2(t)  +  c*3(t)  . 


(32) 

(33) 

(34) 


Consequently,  K^(t)  is  called  the  effective  plane  strain  relaxation  bulk 
modulus. 


If  we  assume  the  average  state  of  stress  to  be 


°22^  =  533^  "  »  (35) 

Oi3(t)  =  a12(t)  =  ^(t)  =  0  *  (36) 

and  retain  the  plane  strain  condition 

en(t)  =  0  ,  (37) 

then  it  is  easily  shown  that 

e22(t)  =  e33(t)  =  e  (t)  ,  (38) 

1  f1-  *  dnfx') 

c(t)  =  y  kT(t-x)  —  —  di  ,  (39) 

o 


where  luj,(t)  |s  the  effective  plane  strain  bulk  creep  compliance.  It  is 
related  to  by 

piL(p)  -  —rj -  .  (40) 

pkT(p) 

Now  assume  the  mixed  state  of  average  stress  and  strain  to  be 

En(t)  *  c(t)  ,  (41) 

o22(t)  =  °33(t)  =  0  >  (42) 

and  all  other  average  shear  stresses  and  strains  vanish.  This  situation 
corresponds  to  uniaxial  straining  in  the  fiber  direction  of  a  cylinder 
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with  free  lateral  surface.  Under  these  circumstances,  the  relation  be¬ 
tween  uniaxial  average  stresses  and  strains  o  (t)  and  e^(t)  is  given 
by 


5n(t)  ' 


*  dtll(T) 

Ea(t_T)  “ -  dT  • 


(43) 


where  Efl(t)  is  the  effective  axial  Young's  relaxation  modulus.  The  in¬ 
version  of  the  relation  (43)  is 


=u(t) 


*  do  (t) 

ea(t"T)  "IT-  dT  ’ 


(44) 


where  e  (t)  is  the  effective  axial  Young's  creep  compliance.  The  func¬ 
tions  o?  E*  and  ea  are  related  in  the  Lapla  -a  transform  domain  by 


pMp) 


peg(p) 


(45) 


For  discussion  concerning  the  Poisson's  effect  associated  with 
uniaxial  stressing,  see  Reference  [3]. 

We  shall  now  proceed  to  discuss  some  explicit  results  obtained 
in  Reference  [3]  for  effective  relaxation  moduli  and  creep  functions  of 
viscoelastic  fiber  reinforced  materials.  It  is  again  recalled  that  all 
results  are  based  on  the  composite  cylinder  assemblage  model. 

It  is  assumed  that  the  fibers  are  isotropic  elastic.  The  fi¬ 
ber  elastic  moduli  are  given  the  subscript  f.  Various  moduli  to  be  used 
are 


\ 


f 


Lam4  modulus , 


Af-KJf 


Shear  modulus , 

Poisson's  ratio, 

Young's  modulus. 

Plane  strain  bulk  modulus. 


It  is  assumed  that  the  matrix  is  isotropic  viscoelastic.  In 
order  to  write  the  matrix  stress-strain  relations,  let  the  stresses  and 
strains  first  be  split  into  isotropic  and  deviatoric  parts.  Thus, 


o 


ij 


+  s 


ij 


■  0  -  1/3  °kk  • 


(46) 
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[ij  ■  c6ij  +  eiJ  '  1  *  1/3  'kk  • 


(47) 


The  most  general  viscoelastic  stress-strain  relations  at  any 
point  in  the  matrix  have  the  form 


o(t) 


(48) 


Sij(t)  *  2 


t  de  (t) 

Gm:t-T>  -  dl  • 

o 


(49) 


E  (t)  =  1/J 


I_(t-T)  dT  , 


m 


(50) 


e  (t)  =  1/2 


ft  ds^x) 

J  (t-x) 


il—  dx 


dx 


(51) 


Here,  K^t)  and  Gm(t)  are  the  matrix  bulk  and  shear  relaxation 
moduli,  respectively,  and  Im(t)  and  Jm(t)  are  the  matrix  bulk  and  shear 
creep  compliances,  respectively.  The  relaxation  moduli  and  creep  com¬ 
pliances  are  connected  through  their  Laplace  transforms,  as  in  previous 
examples,  by  the  relations 


pKm(p) 


pIm(p) 


(52) 


pGm(p) 


pJm<P> 


(53) 


The  matrix  properties  as  expressed  by  the  functions  K^t), 

Gm(t),  Im(t) ,  and  Jm(t)  may  be  regarded  as  experimentally  measured  informa¬ 
tion.  It  is  also  customary  to  approximate  viscoelastic  stress- strain 
relations  by  time  differential  operator  relations  which  may  be  described 
by  spring-dashpot  models.  Although  these  operator  relations  are  a  less 
satisfactory  description  than  the  general  relations  (48)  through  (51), 
they  are  used  for  reasons  of  mathematical  convenience.  It  should  be 
noted  that  (48)  through  (51)  include  spring-dashpot  model  relations  as 
special  cases. 

After  Reference  [3],  we  introduce  for  reasons  of  convenience 
the  following  notation: 


K 

m 


(P)  =  P^P)  » 


(54) 


9 


(55) 


'm(p>  '  pCn(p)  • 


expression : 


The  Laplace  transform  of  k^(t)  is  then  given  by  the  following 


pKT(p)  =  Km(p)  +  3  rm(p)  + 


1  -  C 


h  -  vp)  - 1  r„<p> 


*m(p)  +  3  'm<p) 


-1 

c 


(56) 


where  c  is  the  fiber  volume  fraction. 

Frequently  it  is  permissible  to  assume  that  a  viscoelastic 
material  behaves  elastically  for  isotropic  stress  and  that  the  visco¬ 
elastic  effect  is  confined  to  shear.  Then  <m(p)  in  (56)  is  replaced  by 
an  elastic  bulk  modulus  of  the  matrix,  K|j.  Also,  in  practice,  the  fi¬ 
bers  are  usually  very  much  ctiffer  than  the  matrix.  In  such  cases,  the 
fibers  may  be  assumed  to  be  rigid.  This  last  simplification  when  used  in 
(56)  yields  the  simple  result 

Vp>  -  Vp>  +  i  e„(p)  +  ‘ Vp>  + 1  Vp>'  rhr  •  (57> 

which  may  be  directly  inverted  to  the  time  domain  and  then  becomes 

K*(t)  =  K  (t)  +  ~  G  (t)  +  [K  (t)  +  4  G  (t)  ]  .  (58) 

I  m  Jm  m  dm  1-c 

The  Laplace  transform  of  Gg(t)  is,  according  to  Reference  [3], 
„*  G,(l  +  c)  +  pG  (p)(l  -  c)  „ 

G  (P)  - - : -  G  (p)  .  (59) 

Gf (1  -  c)  +  pGm(p)(l  +  c) 

If  the  fibers  are  assumed  to  be  rigid  relative  to  the  matrix, 
(59)  simplifies  to 


G>  ■  rH  a„(p>  •  <60) 

which  can  be  directly  inverted  to  the  time  domain  to  read 

G*(t)  =  G  (t)  .  (61) 

a  1  -  c  m 

The  Laplace  transform  of  E*(t)  is  given  approximately,  but 
with  a  high  degree  of  accuracy,  by 
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(62) 


E*(p)  =  (1  -  c)  En(p)  +  ^  Ef  , 

which,  when  inverted  to  the  time  domain,  reads 

E*(t)  =  (1  -  c)  E  (t)  +  c  E  H(t)  ,  (63) 

a  m  r 

where  H(t)  is  the  Heaviside  unit  step  function. 

Results  for  the  viscoelastic  Poisson's  effect  and  the  shear 
relaxation  modulus  Gj(t)  are  more  complicated  and  will  not  be  discussed 
here.  For  discussion,  see  Reference  [3]  and  subsequent  developments  in 
this  report. 
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3.  VIBRATIONS  AND  COMPLEX  MODULI  OF  HOMOGENEOUS, 
ANISOTROPIC  VISCOELASTIC  BODIES 


We  shall  begin  with  a  discussion  of  the  concept  of  complex 
moduli  for  anisotropic  and  homogeneous  viscoelastic  materials.  The  dy¬ 
namic  equations  of  such  materials  in  the  absence  of  body  forces  may  be 
written  in  the  following  form: 


,2 

9  u . 

niJ’J  =  p7/ 


9  e.„(x,T) 

Cijk2^t-T)  9t  dT  ’ 


(64) 


(65) 


CU  '2  <“l.J  +  “J,i> 


(66) 


Here  ojj  ,  u^,  and  ejj  are  the  local  stresses,  displacements,  and  small 
strains,  respectively;  p  is  the  density;  and  are  the  relaxation 

moduli.  It  is  seen  that  there  are  6  stresses,  3  displacements,  and  6 
strains  (15  quantities  in  all)  .  Now  there  are  3  equations  (64)  ,  6 
equations  (65)  ,  and  6  equations  (66)  (15  equations  in  all) .  The  stresses 
and  strains  can  be  easily  eliminated  from  (64)  through  (66)  by  substitu¬ 
tion  of  (66)  into  (65)  and  the  resulting  expressions  into  (64).  Thus, 
one  obtains  3  equations  for  the  three  uj.  This  will  not  be  done  here 
since  the  resulting  equations  are  very  complicated  and  will  not  be 
needed  here. 


With  Equations  (64)  through  (66)  there  must  be  associated 
boundary  and  initial  conditions.  A  general  kind  of  such  conditions  for 
a  body  of  volume  V  and  surface  S  is  given  by 


u1(S,t)  = 

T^S.t)  = 

u^x.o)  = 
9ui 

— U,o>  . 

Here,  T^  are  tractions  expressed  b 


o 
u , 

on  S  , 

(67) 

i 

u 

on  ST  , 

(68) 

f  ±  (ii)  , 

(69) 

gi(x)  • 

(70) 

T. 

l 


°ijnj 


(71) 
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and  on  the  surface  S,  n^  are  the  components  of  the  outward  unit  normal. 
The  displacements  are  prescribed  by  (67)  as  functions  of  time  and  surface 
coordinates  on  the  part  Su  of  the  surface.  Similarly,  the  tractions  are 
prescribed  on  the  part  Sp  of  the  surface;  of  course,  Su  +  Sj  =  S.  Con¬ 
ditions  (69)  and  (70)  express  the  fact  that  the  displacements  u-^  and  the 
3u£ 

velocities  are  given  at  initial  time  zero. 

We  do  not  wish  to  look  into  the  dynamic  problem  in  such  gener¬ 
ality.  Instead,  we  consider  the  problem  of  steady-state  vibrations.  In 
this  case  the  displacements  have  the  form 


U.  (x,t)  =  ui(x)  eu0t  .  (72) 

Here,  u^(x)  are  three  space-variable  functions,  \  -  /-I,  and  w  is  the 
circular  frequency.  Insertion  of  (72)  into  (661  yields 


c.j(x)  =  c4  (x)  e 


\  -ot 


where 


Cij  =  2  (ui,j  +  uj,k)  ’ 
If  we  insert  (73)  into  (65),  we  find  that 


o . . (x,t)  =  c,  (x)  iw 


ij 


k 


CijU(t-T)  el“T  dt 


(73) 


(74) 


(75) 


However,  it  should  be  noted  that  a  physical  approximation  is 
involved  in  this  substitution,  for  (65)  is  based  on  the  Boltzmann  su¬ 
perposition  principle  in  time  and  is  therefore  limited  to  gentle  time 
variations.  This  leads  us  to  the  conclusion  that  (75)  is  valid  for 
vibrations  which  are  not  too  rapid.  Since 


a)  =  2nn  , 


(76) 


where  n  is  the  number  of  cycles  per  unit  time,  the  validity  of  (75)  must 
be  limited  to  small  frequency. 


It  is  convenient  to  perform  the  change  of  variable 

i  '  =  t  -  i 

in  (75)  .  Then  the  equation  transforms  into 

-t 


°ij  Qi.t)  =  ck£(x)  e 


iwt 


CiJk£<I,)e"l“T  * 


(77) 


(78) 
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The  time  dependence  on  the  right  side  of  (78)  is  not  a  sinusoidal  vibra¬ 
tion,  for  besides  the  eiwt  factor  the  integral  is  dependent  on  t  because 
of  the  upper  limit  of  the  integral.  The  assumptions  are  now  made  that  t 
is  so  large  that  it  can  h"  given  infinite  value  in  the  upper  limit  of 
the  integral  and  that  thxs  integral  then  converges.  Note  that  e'ut  is 
not  affected  by  making  t  large,  because  it  is  purely  oscillatory  for  all 
times.  Then  we  find  that 


(x,t)  =  ck£(x)  e 


Kilt 


(79) 


We  introduce  the  notation 


IU) 


CijkK(T')  e‘'“T  d,‘  ‘  D1Jkli<,")  • 


(80) 


Then  (79)  can  be  written  in  the  form 


1,(k,t)  =  o  <x)  e1Lot  , 

(81) 

ij  (— '  ■  Dljkt(lu)  ‘kl(2>  • 

(82) 

It  is  seen  that  now  the  stress  cime  variation  is  of  the  nature  of  a 
sinusoidal  oscillation  and  that  this  happens  only  after  a  very  long  time 
has  elapsed.  This  is  physically  reasonable  since  a  steady  state  of 
vibrations  can  be  achieved  only  after  sufficient  time  has  elapsed  for 
transient  effects  to  die  out. 


Because  of  the  formal  resemblance  of  (82)  to  an  elasticity 
stress-strain  relation,  the  quantities  Dijk£  are  called  the  complex  mod¬ 
uli.  They  may  also  be  written  in  the  form 


Dijkn<'“)  ■  Dijki(lo)  4  ,Dijk/(<“)  • 


(83) 


R  I 

where  D  and  D  ,  and  the  real  and  imaginary  parts,  respectively,  of 

ijKX.  ijK£ 

i  j  k£  ’ 

If  (81)  and  (72)  are  substituted  into  (64),  the  exponential 
cancels  out  and  we  find  that 


+  pu 


(84) 


It  is  now  seen  that  the  space-dependent  parts  of  the  stresses,  displace¬ 
ments,  and  strains  enter  into  Equations  (74),  (82),  and  (84).  These  are 
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i 


t 

i 

i 


precisely  the  equations  governing  the  space-dependent  parts  of  these 
quantities  in  the  problem  of  elastic  vibrations,  with  elastic  moduli 
Ci  1  Rx,  replaced  by  complex  moduli  D-j^ j k.£  • 


We  still  have  to  discuss  boundary  and  initial  conditions  for 
the  space-dependent  parts.  For  this  purpose,  (67)  and  (68)  have  to  be 
taken  in  the  oscillatory  form 


u1(S,t)  - 
Tj  (S ,  t)  =  T° 


<S)  e'ut 

on  S  , 
u 

(85) 

(S)  e'”1 

on  ST  . 

(86) 

In  view  of  (71),  (72),  and  (73),  the  boundary  conditions  (85)  and  (86) 
assume  the  form 

u^S)  =  u°  ,  (87) 

T^S)  =  T°  .  (88) 

With  these  boundary  conditions,  the  analogy  with  the  elastic  vibration 
problem  is  complete.  Equations  (74),  (82),  and  (84)  together  with  (87) 
and  (88)  define  a  unique  solution  to  the  problem.  No  initial  conditions 

can  be  prescribed.  This  is  logical,  for  the  assumed  oscillatory  nature 

of  all  variables  prescribes  the  time  variation.  Furthermore,  this  time 
variation  is  steady  state;  it  has,  so  to  speak,  "neither  a  beginning  nor 
an  end",  and  therefore  initial  conditions  cannot  be  imposed. 


In  practice,  boundary  conditions  of  type  (85)  and  (86)  will  be 
given  either  in  the  form 


u .  (S ,  t)  =  u,  (S)  cos  u>t  on  S 


T^S.t 

or  in  the  form 

ui(S,t 

Ti(S,t 

It  is  recalled  that 


u 


T°(S) 

cos 

U)t 

on 

s  , 

T 

u°(S) 

sin 

U)t 

on 

Su  * 

T°(S) 

sin 

(Alt 

on 

ST  • 

U)t  = 

Re(e 

ICOt^ 

> 

n 

4-1 

3 

Im(e 

lUlt^ 

(89) 

(90) 


(91) 

(92) 


(93) 
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Since  the  coefficients  in  (92)  are  complex,  the  solution  for 

and  o^j  will  also  be  complex.  If  the  boundary  conditions  are  of  type 
(89)  and  (90),  the  real  part  of  the  solution  is  taken.  If  the  boundary 
conditions  are  of  type  (91)  and  (92),  the  imaginary  part  of  the  solution 
is  taken. 


In  view  of  the  similarity  to  the  elastic  vibration  problem,  we 
can  now  formulate  an  analogy:  To  solve  a  viscoelastic  vibration  prob¬ 
lem,  take  the  solution,  for  the  corresponding  elastic  vibration  problem. 
In  this  solution  replace  the  elastic  moduli  by  the  complex  moduli.  The 
resulting  solution  has  a  real  part  which  corresponds  to  boundary  condi¬ 
tions  of  type  (89)  and  (90)  and  an  imaginary  part  which  corresponds  to 
boundary  conditions  of  type  (91)  and  (92). 


In  view  of  further  developments,  it  is  important  to  point  out 
an  interesting  relation  between  the  relaxation  moduli  and  the  complex 
moduli.  Let  the  Laplace  transforms  of  the  relaxation  moduli  C^j^^t)  be 
defined,  as  in  Section  2,  by 


^ijkH(p) 

Furthermore,  the  LT  of  (65)  is 

°ij  = 

In  analogy  with  the  definition 
here  the  local  TD  moduli 


■00 


'  o 
by 


P  c 


Cijk£(t)  6  P  dt‘ 

the  convolution  theorem 


(94) 


(95) 


(6)  for  effective  TD  moduli,  we  define 


Bijkt(p)  '  p  Cijkl(p>  ’  <96) 

If  we  now  compare  (96)  in  terms  of  (94)  with  (80) ,  we  note  the  very 
important  relationship 


Dijk«<,“)  '  “ijkt'"*0 


(97) 


That  is:  If  the  TD  moduli  are  known  as  functions  of  the  LT  variable 
then  the  complex  moduli  are  given  by  replacement  of  p  by  iw. 


P> 
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4,  EFFECTIVE  COMPLEX  MODULI  OF 
UN  I AXIALLY  FIBER  REINFORCED  MATERIALS 


The  extension  of  the  theory  which  has  been  developed  in  Sec¬ 
tion  3,  to  composite  bodies  such  as  fiber  reinforced  materials  is  not  a 
simple  matter.  To  try  to  find  a  detailed  solution  of  the  dynamic  equa¬ 
tions  in  fibers  and  matrix  is  an  absolutely  hopeless  task,  and  any  treat¬ 
ment  must  be  limited  to  utilization  of  effective  physical  properties. 

The  questions  which  arise  now  are:  is  it  permissible  to  use  effective 
physical  properties  in  the  equations  developed  in  Section  3,  and  will 
the  results  of  the  calculation  yield  the  macroscopic  dynamic  behavior  of 
viscoelastic  fiber  reinforced  specimens?  The  answer  to  these  questions 
is:  theoretically  no  but  practically  yes,  if  the  frequencies  of  vibra¬ 

tion  are  not  high. 

To  discuss  this  statement,  it  should  be  remembered  that  calcu¬ 
lation  of  static  effective  properties  such  as  effective  relaxation 
moduli  is  based  on  statistically  homogeneous  states  of  stress  and  strain 
[3,4].  Such  states  of  stress  and  strain  occur  in  fiber  reinforced  mate¬ 
rials  in  static  equilibrium  for  certain  boundary  conditions,  which  are 
called  homogeneous  boundary  conditions  [3].  However,  in  dynamic  prob¬ 
lems,  statistically  homogeneous  states  of  stress  and  strain  apparently 
can  not  occur,  because  even  for  homogeneous  bodies  there  does  not  exist 
a  single  case  of  a  dynamic  solution  in  which  the  states  of  stress  or 
strain  do  not  vary  in  space. 

On  the  other  hand,  experience  teaches  us  that  if  continuum 
theory  with  effective  properties  is  applied  to  dynamic  problems,  for 
moderate  frequencies,  the  results  predicted  approximate  quite  well  glob¬ 
al  quantities  such  as  beam  deflections  and  resonant  frequencies.  The 
investigation  of  the  limitations  of  this  approximation  is  extremely 
difficult,  and  to  date  these  limitations  have  not  been  explored.  There¬ 
fore,  it  will  be  assumed  in  the  present  work  that  continuum  theory  with 
effective  properties  is  meaningful  for  fiber  reinforced  composites. 

Sore  justification  for  this  may  be  provided  by  the  following 
argument.  It  may  be  imagined  that  a  fiber  reinforced  specimen  is  di¬ 
vided  into  many  representative  volume  elements  (RVE).  Each  RVE  is  a 
small  part  of  the  specimen,  yet  its  cross  section  contains  many  fibers. 
The  stresses  and  strains  in  such  an  element  may  be  thought  to  be  locally 
statistically  homogeneous.  Then  the  relation  between  the  average  stress 
and  strain  in  any  element  is  given  by  (1)  and  (2),  although,  in  contrast 
to  the  static  case,  the  averages  vary  from  one  RVE  to  another.  The 
equations  of  motion  of  an  element  then  assume  a  form  similar  to  (66) 
through  (66)  with  partial  space  derivatives  replaced  by  increments  of 
quantities  over  space  increments,  i.e.,  sides  of  cubic  RVE.  The  result¬ 
ing  difference  equations  are  approximated  by  the  differential  equations 
(64)  through  (66). 


17 


Adoption  of  the  fundamental  assumptions  discussed  above,  in 
conjunction  with  the  results  given  in  Section  2  for  effective  relaxation 
moduli,  leads  directly  to  calculation  of  the  effective  complex  moduli. 
The  transverse  isotropy  of  a  uniaxially  fiber  reinforced  material,  with 
respect  to  the  x^  axis  in  fiber  direction,  simplifies  (82)  in  the  same 
way  that  (8)  simplifies  to  the  form  (10)  through  (15),  as  a  consequence 
of  the  same  symmetry.  We  can  therefore  write  out  (82)  in  the  following 
form 


*  _ 


*  ^ 


*  ~ 


°11  °11  C11  +  D12  E22  +  °12  E33  » 


° 22  =  °12  £11  +  °22  c22  +  °23  £23  » 


°33  °12  E11  +  °23  E22  +  °22  E33  ’ 


512  =  2  °44  E12 


°13  =  2  °44  e13  » 


(98) 

(99) 

(100) 

(101) 

(102) 


_  *  *  _ 

°23  —  ~  0oa)  c. 


22  23  23 


(103) 


Here,  and  are  the  space-dependent  parts  of  the  local  averages 
a  and  e ,  which  are  of  the  oscillatory  form 


5ij  ■  5ij (i>  c 


uot 


(104) 


tij  =  ^(x)  e  Wt  •  (105) 

The  D*  coefficients  in  (98)  through  (103)  are  the  effective  complex 
moduli.  In  the  DijiC£  notation  appearing  in  (82),  they  have  the  follow¬ 
ing  meaning: 


D11 

=  Dllll  ’ 

(106) 

* 

*  * 

°12 

=  °1122  =  °1133  ’ 

(107) 

* 

*  * 

d22 

=  °2222  =  °3333  ’ 

(108) 

k 

d23 

* 

=  d2233  , 

(109) 
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(110) 


*  *  * 

°44  =  °1212  "  °1313  ' 

Denoting  the  left  sides  of  (106)  through  (110)  by  D^n  and  using  the 
result  (97),  we  have 

D  (iui)  =  B  (ico)  .  (Ill) 

mn  mn 

The  meaning  of  (111)  is  as  follows:  To  find  the  effective  complex 
moduli  D*n,  it  is  necessary  to  replace  the  variable  p  in  the  TD  moduli 
by  the  variable  iw.  The  Bmn(p)  are  known  on  the  basis  of  the  com¬ 
posite  cylinder  assemblage  model  and  have  been  discussed  in  Section  2. 
Thus,  we  have  a  straightforward  method  to  give  expressions  for  the 
effective  complex  moduli. 

It  is  also  quite  obvious  that  we  can  define  physically  mean¬ 
ingful  complex  moduli,  such  as  axial  and  transverse  complex  shear  moduli, 
axial  and  transverse  complex  Young's  moduli,  and  a  plane  strain  complex 
bulk  modulus,  in  analogy  to  similar  moduli  defined  in  Section  2.  We 
shall  give  these  moduli  the  same  notation  as  the  static  relaxation 
moduli  with  an  added  v  above  them.  Thus,  in  analogy  with  (24),  (25), 
and  (34) , 


V  * 

G  (nil) 
a 

"  D*4<lw> 

• 

(112) 

V  * 

* 

★ 

2Gx(iu) 

=  D22(;uj) 

-  D23(iu0  , 

(113) 

v  * 

i k 

* 

3 

CM 

=  D22(i<u) 

+  D23  (  i  co)  . 

( 114) 

The  complex  axial  Young's  modulus  is  defined  for  states  of 
stress  and  strain  analogous  to  (41)  and  (42) .  We  take  a  cylinder  and 
apply  average  uniaxial  oscillatory  stress 

4  \u)t  ,  _v 

all  =  °11  6  ’  (115) 

producing  uniaxial  oscillatory  average  strain 

E11  =  CH  e  »  (H6) 

then 

°U  =  E*(uu)  .  (117) 

The  explicit  expressions  of  the  effective  complex  moduli  of 
the  fiber  reinforced  mateiial  will  involve  the  complex  moduli  of  the 
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fibers  and  the  matrix.  Therefore,  the  meaning  of  these  will  be  briefly 
discussed.  The  fibers  being  assumed  perfectly  elastic  and  isotropic, 
their  complex  moduli  are  simply  their  isotropic  elastic  moduli.  The 
stress-strain  relation  for  fiber  material  is 


’ij  ‘f  Lkk  ij  +  2Gf  cij  • 


If  o  and  e  have  the  oscillatory  forms 


’ij  =  °ij  6 


lut 


(118) 


(119) 


not 


(120) 


then  insertion  of  (119)  and  (120)  into  (118)  gives  simply 


ij 


=  \  i  6 

f  kk  ij 


+  2G , 


(121) 


We  see  that  the  elastic  moduli  Af  and  Gf  are  trivially  the  complex  mod¬ 
uli.  The  physical  meaning  of  this  is  that  in  a  perfectly  elastic  body 
oscillatory  stresses  and  strains  are  in  phase. 


For  the  viscoelastic  matrix  the  situation  is  of  course  differ¬ 
ent.  It  has  been  described  in  detail  for  a  general  viscoelastic  aniso¬ 
tropic  material  in  Section  3,  and  the  present  isotropic  material  is  a 
special  case.  Going  through  the  same  arguments  which  lead  to  (82),  we 
find  for  the  present  isotropic  viscoelastic  matrix  a  relation  of  the 
form 


o.,  =  A  (icj)  e  6  +  2G  (iw)  e  .  (122) 

ij  m  kk  ij  m  ij 

Here,  o^j  and  eij  are  the  time-independent  parts  of  (119)  and  (120). 
^m(ujj)  and  CtoOoj)  are  the  complex  moduli  of  the  matrix,  and,  since  they 
are  functions  of  iw,  the  stresses  and  strains  are  out  of  phase.  The 
complex  moduli  of  the  matrix  must  be  determined  by  a  vibration  experi¬ 
ment,  in  which  a  specimen  is  subjected  to  oscillatory  stresses  and 
strains  with  various  frequencies  w.  Note  that  o^j  and  e^j  are  the 
amplitudes  of  the  stress  and  strain  variations.  The  relation  between 
these  yields  the  complex  moduli. 

In  analogy  to  elasticity,  various  other  related  complex 
moduli  can  be  defined;  e.g., 


E  ( iw) 
m 


V 

G 


m 


+  2G  ) 
m 


V 

A 

m 


V 

+  G 

m 


(123) 
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(124) 


V  V 

K  (iw)  =  X 
m 


m 


2  v 
+  T  G 
3  m 


k  (  \oj) 
m 


V 

+  G 

m 


» 


(125) 


where  (123)  is  the  complex  Young's  modulus,  (124)  is  the  complex  three- 
dimensional  bulk  modulus,  and  (125)  is  the  complex  plane  strain  bulk 
modulus.  Generally,  the  complex  Young's  modulus  is  the  easiest  to  mea¬ 
sure,  and  the  complex  shear  modulus  can  also  be  determined  experimen¬ 
tally.  All  the  other  related  complex  moduli  can  then  be  computed. 

We  now  proceed  to  construct  detailed  expressions  for  some  of 
the  effective  complex  moduli  of  a  uniaxially  fiber  reinforced  material 
by  use  of  the  method  outlined.  We  start  with  the  effective  plane  strain 
bulk  modulus  as  given  by  (50),  and  according  to  what  has  been  said 
before,  the  transform  variable  p  is  replaced  by  iu  everywhere.  We  then 
have 


v  * 

K„(iuj)  =  <  (uj) 
I  m 


4 


r  (ui) 

m 


1  -  c 


4 

k  -  <  (u»)  -  -r  I’  (iw)  <  (iu)  +  cr  r  (uo) 
f  m  3  m  m  3  m 


-1 


c  .(126) 


The  left  side  of  (56)  is  the  effective  TD  plane  strain  bulk  modulus. 
Therefore,  according  to  (96)  and  (97),  the  left  side  of  (126)  is  the 
effective  plane  strain  complex  modulus  k£(u±j).  Also,  ^(p)  and  rm(p)  in 
(56)  are  TD  bulk  and  shear  moduli,  respectively,  of  the  matrix.  There¬ 
fore,  Km(iw)  and  rm(iw)  are  the  complex  bulk  and  shear  moduli,  respec¬ 
tively,  of  the  matrix;  j.e., 

k  (uu)  -  K  (uj)  ,  (127) 

m  m 

r  (iw)  =  G  (\u)  .  (128) 

m  m 


So  (126)  assumes  the  form 


K  (\uj)  =  K  (iw)  +  \  G  (iw) 
1  m  3  m 


_ 1 _ 

K  (iu)  -  4  G  (iw) 
m  3  m 


+ 


1  -  c 


K  (iu>) 

ii. 


V 


G  (uj) 
m 


-1 

c  .(129) 
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When  the  fibers  are  assumed  to  be  rigid  in  comparison  to  the  matrix, 
which  is  a  very  good  approximation  for  a  fiber  reinforced  material,  we 
multiply  (57)  by  p  on  both  sides.  Taking  into  account  (54)  and  (55), 
we  then  have 

K*0W)  =  <  (uo)  +  \  r  (loo)  +  [k  (iu>)  +  4  i  (iw)]  .  (131) 

1  m  3  m  m  3  m  1-c 

Again  the  left  side  is  K^(iw),  and  in  view  of  (127)  and  (128),  we  can 
rewrite  (130)  in  the  form 

K*(xu>)  =  K  (tu)  +  ^  G  (iu)  +  [K  (iu)  +  4  G  (iu>)]  .  (131) 

T  m  dm  m  3m  1-c 


It  is  often  convenient  to  split  a  complex  modulus  into  its 
real  and  imaginary  parts.  Thus,  we  may  write 


v  * 


m 


*R. 


<*< 

m 

3*< 

m 


*1 

^  (co)  , 

(132) 

K^(cu)  , 
m 

(133) 

G*(u>)  . 
m 

(134) 

It  is  seen  that  (131)  separates  very  simply  into  real  and  imaginary 
parts.  Thus, 


k*r  =  KR  +  ~  gr  +  [kr  +  y  gr]  , 

T  m  3m  m  3  m  1  -  c 

K*1  «  K1  +  |  G1  +  [K1  +  ^  G1]  t-5—  . 
T  mjm  m  jml-c 


(135) 

(136) 


On  the  other  hand,  the  separation  of  (129)  into  real  and  imaginary  parts 
is  quite  complicated.  It  is  seen  that  and  will  be  functions  of 

all  of  the  real  and  imaginary  parts  of  ^(iui)  and  £m(iui). 

V  £ 

We  now  turn  to  the  effective  axial  shear  modulus  Ga(iuj).  We 
multiply  (59)  on  both  sides  by  p.  Taking  into  account  (55),  the  follow¬ 
ing  expression  is  obtained: 


P 


V  * 

G  (p) 


g  (l  +  c)  +  r  (p) (i  -  c) 

f  FTC  T-.  /  \ 

G,(l  -  c)  +  T  (p)(l  +  c)  m'P'  • 
r  m 


(137) 


Now  p  is  again  replaced  by  iw.  Then  the  left  side  of  (137)  becomes 

fi*(iuj),  and  in  view  of  (128),  (137)  is  rewritten  as 
a 


22 


(138) 


G  (1  +  c)  +  G  (\us)  (1  -  c) 

f  -  G  (iiii)  . 

m 


v*  f 

Vlu)  „  - - - 

Gr(l  -  c)  +  G  (iu) (1  +  c) 


m 


If  the  fibers  are  assumed  to  be  rigid,  then  on  the  basis  of 


(60), 


V  *  I  +  c  v  v 

G  (ioj)  =  4 -  G  (lu)  • 

a  i  -  c  m 


(139) 


This  simple  expression  separates  easily  into  real  and  imaginary  parts 
Writing 

*R,  ,  .  „*I 


V  ft  ft  W  ft  I 

G  ( i (jj )  =  G  (w)  +  i  G  (lo)  , 
a  a  a 


(140) 


and  using  (134) ,  expression  (139)  can  be  rewritten  as 


X*R  1  +  c  R 

u  =  ■;  u  , 

a  1  -  c  m 


5*1 .  fii  Gi 

a  1  -  c  m 


(141) 

(142) 


The  shear  loss  angle  6  of  the  matrix  is  defined  by 


tan  6  =  — 

"  gr 

m 


(143) 


Similarly,  an  axial  shear  loss  angle 
by 


* 

tan  <5  = 


a 


6^  can  be  defined  for  the  composite 


(144) 


Inserting  (141)  and  (142)  into  (144)  and  comparing  with  (143),  we  obtain 
the  significant  result 

* 

tan  6  =  tan  6  .  (145) 

am 


Thus ,  the  presence  of  rigid  fibers  does  not  affect  the  loss  angle. 

Discussion  of  the  effective  complex  axial  Young's  modulus  is 
in  all  respects  similar.  Multiplication  of  both  sides  of  (62)  by  p 
gives 

p  E  (p)  =  (1  -  c)  p  E  (p)  +  c  E  (146) 

a  ID  L 


23 


which  yields  for  the  complex  modulus 


E  (lo)  =  (1  -  c)  E  (uo)  +  c  E  . 
a  m  f 

Separation  into  real  and  imaginary  parts  is  straightforward;  i.« 


*R 


(1  -  c)  E  +  c  E  , 
m  t 


*1  I 

E  =  (1  -  c)  E  . 
ci  m 


(147) 

(148) 

(149) 
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5.  LOr :ITUDINAL  VIBRATIONS  OF  VISCOELASTIC 
FIBER  REINFORCED  ROD 


We  consider  now  a  typical  vibration  problem  of  a  rod  which  is 
built  in  at  one  end  and  is  subjected  to  sinusoidal  time-variable  force 
at  the  other  end  (Figure  2).  The  fibers  are  in  the  direction  of  the 
axis  of  the  rod. 


The  appropriate  end  conditions  are 
u  (o,t)  =  0  , 


o  (l.t)  = 


=  o  e 
x=^  ° 


not 


(15C) 

(151) 


where  u  (o,L)  is  the  longitudinal  displacement,  w  is  the  frequency,  and 

P 

(152) 


no  A 


where  A  is  the  cross-sectional  area.  In  accordance  with  the  general 
discussion  in  Sections  3  and  4,  we  consider  first  the  associated  elastic 
problem.  The  governing  differential  equation  is  the  one-dimensional 
wave  equation 


32u 

9x2 


1  9  u 

2  o  2  * 
c  9 1 


(153) 


where 


(154) 


in  which  is  the  effective  axial  Young's  modulus,  and  p  is  the  average 
density  of  the  material.  We  set 


u  (x,t)  =  U(x)  ela)t  . 


(155) 


Substitution  of  (155)  into  (153)  gives  the  ordinary  differen¬ 
tial  equation 


=  0  . 


(156) 


Substitution  of  (155)  into  (150)  and  (151)  gives  the  end  conditions 

U(o)  =  0  ,  (157) 
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t 

Figure  2.  Longitudinal  Vibrations  of  Viscoelastic  Fiber 
Reinforced  Rod. 


dU(fc)  fo 
dx  * 


The  solution  of  (156)  with  end  conditions  (157)  and  (158)  is 


(158) 


o  c  sin  I 
o  \ 

„  c  / 

*  / 

1  U)i  \ 

wE  cos 

— 

a  \ 

.  c  1 

The  resonant  frequencies  are  found  by  letting 

cos  (^r) =  0 


(159) 


(160) 


so  that  U  becomes  unbounded.  Equation  (160)  in  conjunction  with  (154) 
yields 


w 

n 


(2n  +  1)  tt 
2S 


(161) 


and  the  first  resonant  frequency  is  accordingly 


U) 

o 


IT 

U 


(162) 


According  to  the  theory  in  Sections  3  and  4,  the  modulus  E*  in 

V  * 

(159)  must  now  be  replaced  by  the  effective  complex  modulus  Ea(iuj)  as 
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given  by  (147),  (148),  and  (149).  It  is  convenient  to  represent  ^*(16) 
in  the  usual  form 


v  *  ,  v  *  i  (5 

E  (uj)  =  E  e 
a  1  a  1 


(163) 


where 


i e* i .  ,/  (e!r>2  +  (e!1)2  , 


a 


a 


(164) 


tan  6 


.  S'* 

*  a 


*R  ' 


(165) 


In  view  of  (148)  and  (149), 


i v*  i 

E 

'  a 1 


/(VmEm  +  VfEf)2  +  (VmEm)2  * 
■J  mm  f  f  mm 


(166) 


tan  6 


v  E1 
m  m 


R  * 

v  E  +  v  E 
mm  f  f 


(167) 


where  vm  and  Vf  are  the  volume  fractions  ot  fibers  and  matrix,  respec¬ 
tively. 


A  typical  example  of  a  viscoelastic  fiber  reinforced  material 
is  an  epoxy  reinforced  by  glass  or  boron  fibers.  For  such  materials  it 
is  generally  true  that 


—■  >  «  <  0.1  .  (168) 


The  Ej,,  in  the  left  part  of  (168)  is  the  static  Young's  mod¬ 
ulus,  which  may  be  interpreted  as  E^(io)  at  zero  frequency.  Experiments 
described  in  Reference  [5]  have  shown  that  for  an  epoxy,  E®  is  a  very 
slowly  increasing  function  of  u,  it  being  necessary  to  increase  the 
frequency  from  10  to  10^  in  order  to  produce  a  50%  increase  in  E$. 
Therefore,  the  relative  magnitude  of  Ef  to  E^  may  be  taken  as  about  the 
same  as  that  of  Ef  to  Ejj. 

In  view  of  (168)  and  the  subsequent  remarks,  we  can  approx¬ 
imate  (166)  with  great  accuracy  by 


27 


(169) 


I  V*  R  * 

E  =  v  E  +  vr  E  =  E  (u)  . 

1  a 1  mm  f  f  a 

The  expression  (169)  is  the  effective  elastic  axial  Young's 
modulus  of  a  uniaxially  fiber  reinforced  material  with  the  matrix 
Young's  modulus  Em  replaced  by  E^(w) .  In  the  following  it  will  be 
denoted  by  Eg(oj)  .  It  should  be  noted  that  because  of  (168)  and  the 
subsequent  remarks,  E^(w)  is  very  nearly  equal  to  the  static  effective 
axial  Young's  modulus. 


It  also  follows  from  (167)  and  (168)  that 

* 

tan  6  <<  1  ; 


(170) 


therefore 


*  * 
tan  6  -  & 


(171) 


In  view  of  (163)  and  (169)  , 


J 


E*(  Hu)  .  (eV'V6  11 


(172) 


Define  c  by 


c  - 


E  ( i  w) 
a 


(173) 


From  (154)  and  (172)  , 


c  -  c  e 


*  , 

i(5  n 


Now  replacement  of  E*  by  Eg(iu))  in  (159)  leads  to 


(174) 


U(x)  = 


o  c 
0 


sin 


Of) 


X*,  X  /Uil\  * 

wE^tuu)  cos  I  —  J 

We  use  the  approximations  (172)  and  (174)  in  (175)  to  obtain 


(175) 


o  c  sin 

,  /  N  o  -16  2 

u(x)  =  —I -  e 


E  (oo) 
a 


^  wx  e-1  *5  /2^ 

I  w£  -16  /2\ 
cos  e  ) 


(176) 
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where  u(x)  is  the  space-dependent  part  of  u(x,t)  in  the  viscoelastic 
vibration  problem. 


By  Euler's  formula, 

*  *  * 

-16  /  2  6  .6  /  -1 7  7  \ 

e  =  cos  —  -  \  sin  —  ,  (1/7) 

and  since  by  (170)  and  (171)  6*  is  a  very  small  angle, 

cos  6*/2  -  1  ,  (178) 

sin  6*/2  -  6 * / 2  .  (179) 

(180) 
(181) 

/  oj l  -16  /l\  (  uil  iu£.6  ^ 

cos  (  —  e  )  ~  ccs  y  —  -  i  '2c  j  .  (182) 

Using  the  usual  expression  for  the  sin  and  cos  of  a  complex 
number,  we  obtain 

sin(~  e_l6  /2)  “  sin  [~]  cosh  -  i  cos  sinh  [ii|^,(183) 

cos  l^~  e~x6  /2)  =  cos  cosh  +  i  sin  sinh  (i^'),(l84) 

where  cosh  and  sinh  are  the  usual  hyperbolic  functions.  It  would  now  be 

Co  (5 

very  desirable  to  show  that  is  a  very  small  expression,  so  that  we 

could  use  small  argument  approximations  for  the  hyperbolic  functions. 
Consider  a  typical  glass  fiber  reinforced  epoxy.  Typical  material  prop¬ 
erties  are 

E*  -  0.5  x  106  psi(*) 
m 

E^  =  10  x  10^  psi. 

(*)As  stated  before  is  really  a  very  slowly  increasing  function  of 
frequency.  The  present  constant  value  is  merely  an  order  of 
magnitude. 


Thus , 


-t  c5*/2 


1-16/2  , 


(  cjx  -i  5  /2\  f  wx  (jjxd  \ 

,ln\~e  ]  =  .ln^—  -  i—  ], 
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We  choose  a  large  matrix  loss  angle,  tan  6  =  0.1. 

m 

E1  =  0.05  x  106  psi  . 
m 


Then 


We  take  as  volume  fractions 


v  =  0.  A 
m 


=  0.6 


Then  from  (148)  and  (167), 


E  =  6.2  x  10^  psi 


tan  6  =  0.005  . 

Take  as  typical  density 

p  =  3.0  (relative  to  water). 

Then  it  follows  that 

c  «=  12,770  tt/sec. 

Let  the  length  of  the  rod  be 
£  -  10.0  ft. 

Then 

L  6  9.  _  e  -I  n“6 

— r —  -  2.5  x  10  to  . 

2c 


If 


* 

ui  i 
2z  ' 


<  0.05  , 


we  have  with  great  accur  acy 


sinh 


(1 r) 


*  „  * 

o)6  £ 


2c  » 


cosh 


Tnis  corresponds  to 


o)  <  20,000 


(185) 

(186) 
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Since 


u  =  27m  , 


the  frequency  n  satisfies 


n  <  3180  cycles/sec. 

This  is  a  high  frequency  for  which  the  present  theory  may  not 
be  valid  anymore.  Thus,  it  is  seen  that  (185)  and  (186)  may  be  safely 
assumed. 

Since  o  <  x  <  l,  we  have  with  better  accuracy  than  (185)  and  (186) 


sinh  ( 

t  *  \ 

o)6  x^ 
v  2c  1 

* 

1  0)6  x 

1  "  2c  ’ 

(187) 

cosh  { 

1  oi6  xS 

\  2c  J 

1=1- 

(188) 

Introducing  (185)  through  (188)  into  (183)  and  (184)  ,  then  introducing 
the  resulting  approximations  and  (186)  into  (176),  we  find  the  first 


order  of  magnitude  in 


u)6*x  o)6*£ 

and 


2c 


2c 


u(x) 


(1  -  i  6*/2) 


Since 


u (x , t)  =  u(x)  e1^  , 
we  have  from  (189)  and  (190) 


(189) 


(190) 


o  c  sin 
u  (x ,  t)  =  ~~T - 


*  /  \ 

03  E  cos  — 
a  '  c  ' 


*  A 

[cos  oit  +  6  / 2  sin  oot  +  i(sin  wt  -  6  /2  cos  wt)] 

(191) 


Now  if  the  forcing  stress  (151)  has  a  cosine  variation,  i.e.. 


o(£,t)  =  o  cos  oj t  =  Re  (o  e1ut)  , 
o  o 


(192) 


then  the  displacement  is  the  real  part  of  (191);  thus. 
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(193) 


u£  (x,t) 


sin 

1  01  x  \ 

0  c 

l  —  I 

0 

\  c  l 

1 k 

wE  (w) 

cos  1 

ui£  ^ 

a 

k  c  ; 

(cos  tot  +  6  /2  sin  wt)  . 


If 


o(e,t)  =o  sin  ut  ■  I  (o  ,  (194) 

o  mo 

then  the  displacement  is  the  imaginary  part  of  (191).  Thus, 

o  c  sin  \  * 

u  (x,t)  =  — - 7 — 7—  (sin  OJ t  —  6  / 2  cos  wt)  .  (195) 

ojE  (w)  cos  —  1 
a  V  c  ' 

From  the  elastic  analysis  at  the  beginning  of  this  section, 
it  is  seen  that  the  elastic  solutions  corresponding  to  (192),  (193), 
and  (194),  (195)  are,  respectively, 


u  (x,t)  = 


0  c 
o 


■Ht> 

e*  M  -(7) 


COS  U)t 


(196) 


Ug  (x,t)  = 


0  C 

o 


wEa  (w)  cos[~) 


sin  wt  . 


(197) 


It  is  thus  seen  that  in  the  viscoelastic  vibration  problem , 
there  are  out-of-phase  components  which  are  proportional  to  6*. 

It  is  instructive  to  represent  (193)  and  (195)  in  different 
forms.  It  is  seen  that  these  expressions  may  be  rewritten  as 


(x,t) 


u  (x,t) 
s 


*2 

0  c  1  +  {  74 
°  ¥ 


wE  (10) 
a 


o  c 
o  N 


*2 

1  +  6/4 


wE  (to) 
a 


cos  (wt  -  O  , 


sin  (tot  -  $)  , 


(198) 


(199) 
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where 


* 

tan  <p  =  ~  .  (200) 

Comparison  of  (198)  and  (199)  with  the  elastic  results  (196) 
and  (197)  shows  that  the  amplitude  of  the  viscoelastic  vibrations  is 

increased  by  the  factor  /l  + 5*^/4  and  that  there  is  a  phase  lag  whose 
magnitude  is  given  by  (200).  Since,  as  has  been  shown  before,  6*  is 
very  small  for  a  typical  fiber  reinforced  material,  the  amplitudes  are 
practically  the  same  as  in  elastic  vibrations,  and  also  the  phase  lag 
can  be  safely  ignored.  It  is  also  seen  that  because  of  the  occurrence 
of  cos  (gJ/c)  in  the  denominators  in  (198)  and  (199),  the  resonant 
frequencies  are  the  same  as  for  the  elastic  vibrations  and  are  thus 
given  by  (161)  and  (162). 

The  very  important  conclusion  is  that  for  all  practical  pur¬ 
poses j  the  viscoelastic  effect  can  be  ignored  in  the  problem  of  longi¬ 
tudinal  vibrations  of  fiber  reinforced  rods  made  of  typical  fiber 
reinforced  materials . 

It  should  be  noted  that  this  phenomenon  is  due  to  the  fact 
that  the  fibers  are  very  much  stiffer  than  the  matrix  and  that  the 
matrix  loss  angle  is  small.  In  cases  where  these  conditions  an  not 
fulfilled,  the  foregoing  conclusions  may  become  invalid,  and  it  may  be 
necessarv  to  carry  through  the  analysis  without  the  approximations 
which  wer  introduced  here.  However,  such  cases  do  not  seem  to  be  of 
interest  for  fiber  reinforced  materials. 
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6.  TORSIONAL  VIBRATIONS  OF  VISCOELASTIC 
FIBER  reinforced  cylinder 


As  a  second  typical  problem,  we  consider  a  fiber  reinforced 
cylinder  which  is  built  in  at  one  end  and  is  subjected  to  a  sinusoidal 
torque,  in  time,  at  the  other  end.  Again  we  consider  first  the  associ¬ 
ated  elastic  problem. 


The  basic  variable  is  the  angle  of  twist  0(x,t).  The  govern¬ 
ing  differential  equation  for  0  is 


where 


2 

1  3  6 


2 

c 


* 


G  C 
a 

P  I  * 


(201) 


(202) 


(see,  e.g.,  Reference  [6]).  G*  is  the  effective  axial  shear  modulus, 
and  C  is  a  number  which  depends  upon  the  geometry  of  the  cross  section. 
In  terms  of  the  Saint  Venant  torsion  function  cj>  (y ,  z) ,  C  is  given  by 


C 


(y2  +  z2  +  y  -  z  |£)  dydz  , 
(A) 


the  integral  being  taken  over  the  cross  section,  which  is  referred  to  a 
y,z  coordinate  system  with  origin  at  the  centroid  (see,  e.g.,  [7]). 

Note  that  G*  C  is  the  torsional  rigidity  of  the  cylinder. 

a 

Furthermore,  p  is  the  density  of  the  material  and  I  is  the 
polar  moment  of  inertia  of  the  section. 

The  torque  M(x,t)  is  given  by 


M(x, t) 


The  end  conditions  are 


e(o,t)  =  o  , 


(203) 


(204) 


M(£,t) 


x=£ 


iwt 


e 


(205) 
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In  (205),  M0  is  the  amplitude  of  the  sinusoidal  forcing  torque  and  to  is 
the  frequency. 


We  seek  a  solution  of  the  form 

0(x,t)  =J(x)  ela)t  .  (206) 

Introduction  of  (206)  into  (201),  (204),  and  (205)  yields  the  ordinary 
differential  equation 

,2  2 

^  ^  =  0  (207) 

dx  c 


and  the  end  conditions 

■£(o)  *=  0  , 


d£ 

dx 


x*=£ 


M 

o 

* 

G  C 
a 


The  solution  of  (207),  subject  to  (208)  and  (209), 


is 


(208) 

(209) 


£  (x) 


therefore,  the  complete  solution  is 


0 (x,t) 


C  M 

0 

sin  1 

* 

r  w£  \ 

uj  G  C 

a 

cos  ' 

sin  | 

[— ) 

(  c  /  lU)t 

cos  1 

■ 

The  condition  of  resonance  is  obtained  from 

■  0  . 


then 


cos 


to£  2n  +  1 


(210) 


(211) 


(212) 


c  2 


(213) 


and  from  (202)  and  (213), 
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2n  +  1  Ti 


* 

G  C 
a 


(214) 


iii 

n 


2 


■J  o  I  ’ 


with  first  resonant  frequency 


Li) 

o 


(215) 


We  now  consider  the  case  when  the  nutrix  is  viscoelastic. 

Then  the  space-dependent  part  of  the  angle  0(x,t)  is  denoted  by  0(x)  and 
is  found  by  replacement  of  G*  in  (210)  by  the  effective  complex  axial 
shear  modulus  Qa.  An  expression  for  G  has  been  given  before  for  elas¬ 
tic  fibers  and  viscoelastic  matrix  (138).  We  shall  here  use  the  simpli¬ 
fied  form  (139)  for  rigid  fibers  in  order  to  simplify  the  analysis. 

Since  in  a  typical  fiber  reinforced  material  with  epoxy  matrix  the  ratio 
of  fiber  to  matrix  stiffness  is  about  25,  at  least,  this  is  in  general 
an  excellent  approximation.  In  order  to  avoid  confusion  with  propaga¬ 
tion  velocity  c,  we  shall  denote  the  fiber  volume  fraction  by  Vf.  Then 
we  have 


G  (iu) 
a 


1  +  v 


1 

1  -  m 


f  v 

c  (no)  , 


(216) 


G  (no) 
m 


GR  +  1  G1  , 
m  m 


(217) 


*  m 

tan  6  =  tan  4  =  tan  4  =  — r-  . 

"  CR 

m 


(218) 


Expression  (217)  can  be  rewritten  with  the  help  of  (218)  in 


the  form 


£  (h)  =  G  el*  ,  (219) 

m 

where 

G  =  ,'( GR)2  +  (G1)2  =  /?  W tan2  6  .  (220) 

m  m  J  m 

Hence,  from  (216)  to  (220), 
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(221) 


v*  +  Vf  R  I 

C  (  ioj)  =  - - L  GK  J 

a  l  -  v,.  m 


.  ,  2X  1 6 
1  +  tan  o  e 


We  recall  that  the  elastic  effective  axial  shear  modulus  Ga  is 


given,  in  the  case  of  the  rigid  fibers,  by 


*  1  +  vf 

G  =  - - -  G 

a  1  -  v  m 


(222) 


(see  References  [1  and  3]), where  Gm  is  the  elastic  shear  modulus  of  the 
matrix . 

We  adopt  the  notation 


1  +  Vf  r  * 

"j  ~  G  (w)  =  G  (w)  , 
1  -  v ,  m  a 


(223) 


in  an  analogy  to  the  notation  used  in  (169).  Then  from  (221)  and  (223), 


V  *  £ 

G  (lu)  =  G  ( ou ) 
a  a 


„  2a  i6 

tan  6  e 


(224) 


V  £ 

Denoting  c  as  a  function  ot  G  (iw)  by  c,  we  have  from  (202)  and  (224), 


,  G  (w)  C 
a 

C  P  I 


J 1  +  tan^<5  e^ 


(225) 


Therefore,  from  (202)  and  (225), 


/ •.  .  2  j. .  1/4  1 6/2 

c  =  c(l  +  tan  6)  e 


(226) 


Replacing  c  by  (226)  and  G  by  (224)  in  (210),  we  obtain 


c  M 


-i  6/2 


0  (x)  = 


.  rw  .  2  .-1/4  —  t 6 / 2  . 

sin[—  (1  +  tan  6)  e  x] 
u  _ c _ _ C _ 

w  G  (u>)  C  (1  +  tan^  6)^^  cos[—  (1  +  tan^6) 


(227) 


Since  for  the  usual  viscoelastic  materials  in  general 

tan  6  <  0 . 1  , 

2 

it  is  seen  that  tan  6  can  be  safely  ignored  in  comparison  to  1.  Thus, 
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c  M 


0(x)  = 


ui  G  (co)  C 
a 


.to  -uS/2  . 

,-,6/2  Sl"  (c  6  *> 

cos  (—  e  1 
c 


(228) 


There  is  no  difficulty  in  finding  the  real  and  imaginary  parts 
of  (228).  The  results,  however,  are  very  clumsy,  and  we  shall  therefore 
continue  the  calculation  for  x  =  £ . 

We  have  from  (228)  , 


c  M 


6(0 


to  G  (O  C 

a 


(cos  6/2  -  i  sin  6/2) 


sin[—  (cos  6/2  -  \  sin  6/2)] 


cos[—  (cos  6/2  -  i  sin  6/2)] 
c 


(229) 


We  use  the  notation 


to£.  X  /o 

—  cos  6/2  =  a 
c 

—  sin  6/2  =  6 
c 


(a) 

(b) 


(230) 


Using  the  usual  expressions  for  the  sine  and  cosine  of  complex  quanti¬ 
ties  (see,  e.g.,  Reference  [8]),  we  find 


c  M 


0(0  = 


u)  G  (to)  C 

a 


/  -  ,  r  sin(2a)  -  \  sinh(2B) 

(cos  o/2  -  i  sin  6/2)  - .  V - IVToT  •  (231) 

cos(2a)  +  cosh(28) 


Thus,  the  angle  of  twist  at  x  =  £  is  given  by 


e(£,t)  =  6(0  e 


not  _  C  ^o  sin(2a)  -  \  sinh(2B)  t  (tot-5/2) 

*.  .  coj(2a)  +  cosh(26) 
toC  (to)C 

a 


(23?) 


Tf  (205)  has  the  form 


M(l,t)  =  cos  tot  , 


(233) 


then 


0  ( £  ,  t)  =  Re 


0  ( g.)  e 


itot 


(234) 
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If  (205)  has  the  form 


M(£,t)  =  sin  uit  , 


(235) 


then 


0  U,  t)  =  Im 
s 


r_ 


e(£)  e 


lU)t 


(236) 


It  follows  from  (234)  and  (232)  that 


c  M 


ecU,t)  = 


uG  (w)C 
a 


o  sin(2a)  cos  (cot  -  6/2)  +  sinh(2B)  sin(cit  -  6/2) 

/  n  \  i  t  /On  \  '  \  ^  J  '  / 


cos(2a)  +  cosh(2B) 


From  (236)  and  (232), 


c  M 


„  s  o  sin(2a)  sin  (wt  -  5/2)  -  sinh(2B)  cos(wt  -  6/2) 

°sU,U  ‘  T*7~77  cos(2a)  +  cosh(2B)  * 


uiG  (w)C 
a 


It  is  convenient  to  represent  (237)  and  (238)  in  different  forms.  It 
is  easily  seen  that  (237)  can  be  rearranged  as 


c  M 


ec(2,t)  = 


7 


sin  (2a)  4-  sinh  (26)  ,  .  r  N 

^  cos  (wt  "  »  (239) 


*,  cos(2a)  +  cosh(2B) 

a 


where 


tan  ip 


sinh (26) 
sin(2a) 


(240) 


It  is  thus  seen  that  (239)  is  a  cosine  vibration  with  a  phase  leg  rela¬ 
tive  to  moment  input  (233).  When  (234)  is  written  in  the  form 


©c ( £ , t )  =  Amp  [8c]  cos(wt  -  w )  , 


(241) 


it  is  seen  that  the  amplitude  is  given  by 


c  M 


Amp  [ 0  ^  ] 


o  J~ 


sin^(2a)  +  sinh^(2g) 


*  .  cos (2a)  +  cosh(2B) 

wG  (ijj)C 
a 


(242) 


and  that  the  phase  lag  is  given  by 
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(243) 


'* 


tan 


sinh(2t: ) 
sin(2a) 


+  6/2  . 


Similarly,  (238)  can  be  written  in  the  form 


esU,t) 


c  M 

o 


* 

toG  (uj)C 

a 


/sin^(2ot)  +  sinh^(26) 
cos(2a)  +  cosh(28) 


sin(iut  -  6/2  -  i p)  ,  (244) 


where  ip  is  again  given  by  (240).  Thus,  we  can  also  here  use  the 
representation 

Gs(H,t)  =  Amp  [0  I  sin(o)t  -  , 

where  Amp  [6  ]  is  given  by  (242)  and  ¥  is  given  by  (243). 

Figure  3  shows  a  plot  of  Amp  [6C]  as  a  function  of  input  fre¬ 
quency  for  a  viscoelastic  fiber  reinforced  circular  cylinder.  Also 
shown  is  the  amplitude  variation  with  w  for  the  elastic  vibrations,  with 

matrix  elastic  shear  modulus. 

used  are  as  follows: 

5.0  ft.  length  of  cylinder 

4.0  in.  diameter  of  section 

3.0  density  relative  to  water 

GR(0)  [1  +  t  log, „  ui],  1  <  u  <  104  sec"1 
m  h  1U  -  - 

0.5  x  106  psi 

0.1 

0.4  v^  =  0.6. 

the  neighborhood  of  the  first  elastic  resonant  fre¬ 
quency,  the  viscoelastic  amplitude  has  a  peak  which  is  about  five  times 
that  of  the  initial  amplitude.  The  next  peak  which  comes  behind  the 
second  elastic  resonant  frequency  is  already  very  damped  out.  The  damp¬ 
ing  increases  further  in  the  neighborhood  of  higher  elastic  resonant 
frequencies . 

It  is  thus  seen  that  the  viscoelastic  nature  of  the  matrix  has 
a  very  beneficial  effect  on  torsional  vibrations. 


G^(0)  taken  as  the 

The  data 

l  = 
d  = 

P  *= 

GR  (w)  = 
m 

CR  (0)  = 

m 

tan  6  = 

m 

v  = 
m 

It  is  seen  that  in 
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Figure  3.  Amplitude  of  Angle  of  Twist. 


7.  TRANSVERSE  VIBRATIONS  OF  UNIAXIALLY  FIBER  REINFORCED  BEAM 


We  now  consider  the  problem  of  the  forced  vibrations  of  a  beam 
of  constant  cross  section.  Let  x  be  measured  along  the  beam  axis,  which 
is  also  the  direction  of  fiber  reinforcement.  The  equation  of  motion  of 
an  elastic  fiber  reinforced  beam  is 


,4  i 

3  w  3  w  1  ,  .  . 


(246) 


2 

a 


* 


E  I 
a 

PA  * 


(247) 


where  w  is  the  transverse  deflection,  p  is  the  average  density,  p  is  the 
load  per  unit  length  on  the  beam  as  a  function  of  position  and  time,  I 
is  the  moment  of  inertia  ly,  and  A  is  the  cross-sectional  area. 

We  consider  the  important  case  of  a  simply  supported  beam 
which  is  loaded  at  x  =  £,  by  the  concentrated  force  P0elk)  ,  Figure  4. 

For  this  case  the  boundary  conditions  are 

2  2 

w(o,t)  =  (o,t)  =  w(£,t)  =  (£,t)  «=  0  ,  (248) 

3x  3x 


and  the  loading  function  p(x,t)  assumes  the  form 

p(x,t)  =  Pq6  (x-O  elWt  ,  (249) 


where  6  is  the  delta  function. 

The  solution  to  the  problem  described  by  (246)  through  (249) 
is  well  known  and  is  given  by 


E ,  . 

w  (x,t) 


2P  e 
o 

* 


UiJt 


E  H 
a 


oo 


l 

n=l 


sin(a  O  sin(a  x) 
n  n 


where 


a 

n 


■~TT 

c 


(250) 


(251) 
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Figure  4.  Lateral  Vibrations  of  Viscoelastic  Fiber 
Reinforced  Beam. 


* 

9  /  E  1 

2  4  a 

ui  ■=  ot  — : —  . 

n  n  pA 


We  rewrite  the  solution  (250)  in  the  form 


where 


w  (x,t)  =  elwt  l  S  (x)  . 

n-1  n 


„  2P  sin(a  O  sin(a  x) 

t ,  x  o  n  n 

a  H  7) 

'  n  ' 


In  view  of  (252),  Equation  (254)  can  be  written  as 


v  2P  sin(a  £)  sin(a  x) 

WE(X)  _  — °  - D _ P. . 

nW  £  *  4  2 


E  la  -  ui  pA 
a  n 


The  complex  viscoelastic  solution  then  has  tre  form 


VF  /  .  1U)C  r  _  .  . 

w  (x)  =  e  I  w  (x)  , 


>ut 


n=l 


(252) 


(253) 


(254) 


(255) 


(256) 
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where 


0l  sin(a  O  sin(a  x) 
-  /  v  tt  n  n 

Wn(x)  "  £  v*  a  2 

E  (tus)Iot  -  w  p 
a  n 


In  view  of  (147)  through  (149), 


v*  „  *r  I 

E  (ioj)  -  v  E  +  v.E  +  iv  E 
a  mm  f  f  mm 


We  rewrite  (258)  in  the  form 


(257) 


(258) 


E  ( i to)  -  E  (to)  (1  +  i  tan  6*)  ,  (259) 

a.  3 

where  e£  is  given  by  (169)  and  also  by  (148),  and  tan  6*  is  given  by 
(167).  Substitution  of  (259)  into  (257)  and  rr arrangement  yields 


w.00 


*  4  9  *  4  * 

2P  E  (u))Io  -  uj  p  -  lE  (o))Ia  tan  6 
_ o  a _ n _ _ _ a  n _ 

*  (E  (u)Io4  -  (D2p)2  +  (E  (ca)Ia^  tan  6*) 2 
an  an 


sin  (a  O  sin  (a  x)  . 
n  n 

(260) 


We  shall  adopt  the  notation 


-2 

to 

n 


4  Ea(u)I 
“n  pA 


(261) 


which  is  analogous  to  (252)  .  Note  that  the  right  side  of  (261)  is  a 
function  of  u  because  of  the  presence  of  E*((j).  However,  because  cf  the 
very  slow  variation  of  E*(cj)  with  to  (see  Section  5),  (261)  is  numerical¬ 
ly  quite  close  to  the  value  of  (252).  With  the  notation  (261),  Equation 
(260)  can  be  rewritten  in  the  form 
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Ea(w)I£ 


.  U)  „  * 

1 - 2  ~  1  tan  ” 
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rf  2\ 

2  „  n 

/  ,  0) 

2  * 

1 

LV  “n/ 

+  tan  6 

sin(a  x)  sin(a  O 
n  n 


It  follows  that 


(262) 
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•  sin  (a  x)  sin  (a  £)  .  (263) 

n  n 

This  completely  defines  the  dynamic  viscoelastic  solution  (256).  Thus, 
if  the  applied  force  is  P0  cos  w'_. 


If  the  applied  force  is  Pq  sin  tut,  we  have 
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00  tan  6  sin(a  O  sin(a  x) 
B  -  T  - r-S - S— 


n-1  4 

a 
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i.i  ?  * 

fl - j  +  tan 

u 

n/ 


Equations  (264)  and  (265)  may  be  written  in  the  form 

cos  (tot  -  v?  )  , 

sin  (tot  -  H )  , 


VE 

wc  (x,t)  -  Amp 

L  c  j 

VE 

w  (x,t)  •  Amp 
s 

'  VE 
w 

s 

(267) 


(268) 


(269) 


where 


Amp 


«  Amp 


2P 

o 

E*(to)I£ 

a 


+  b2]1/2 


(270) 


and 

*  -  tan-1  7  .  (271) 

A 

Since  6*  is  a  very  small  angle,  it  is  reasonable  to  neglect 
tan2  6*  in  the  denominator  of  (271),  unless  u>  is  almost  equal  to  ton. 
Barring  this  eventuality,  if  one  remembers  that  Eg  (to)  is  a  very  slowly 
Increasing  function  of  to,  it  appears  that  the  amplitude  is  very  nearly 
that  of  elastic  vibrations  with  Young's  modulus  Eg(0).  Also,  if  io  is 
not  close  to  wn,  the  phase  lag  defined  by  (271)  is  very  small  and  can  be 
neglected. 


It  is  seen  that  the  denominator  in  (270)  cannot  vanish,  and 
thus  there  la  no  resonance  in  theory.  Yet  the  denominator  can  still 
become  very  mall.  Consider  the  case 

u>  ■  £3^  •  (272) 

Since  <Ln  as  defined  by  (261)  involves  E*(uj),  it  is  seen  that  (272)  is  an 
equation  for  ui  which  must  be  solved.  Let  a  solution  of  (272)  be  denoted 
by  tog.  Then  (266)  and  (267)  assume  the  form 
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(b) 


(273) 


It  is  seen  that  for  small  k,  i.e.,  the  first  solutions  of  (272),  the 
last  term  In  (273b)  will  become  very  large*,  thus,  we  have  a  pseudo¬ 
resonance  condition. 

It  is  also  seen  that  when  (272)  is  fulfilled,  Equation  (267) 

gives 

tan  if  -*■  “,  if  «  -|  ,  (27*) 

which  implies  a  phase  reversal  of  the  vibration,  since  then 

cos  (uit  -  if  )  ■  sin  uit 

sin  (uit  -  if )  -  -  cos  uit  .  (275) 

In  order  to  obtain  an  idea  of  the  comparison  of  elastic  and 
viscoelastic  vibrations,  the  deflection  at  the  center  of  a  beam  loaded 
at  its  center  has  been  calculated  numerically.  In  that  case, 

C  -  x  -  ~  ;  (276) 


therefore, 


8in(a  O  sin(a  x)  -  sin 
n  n 


n  odd 
n  even  . 


(277) 


Thus  for  that  case, 
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(279) 


The  amplitude  (279)  has  been  calculated  numerically  as  a  func¬ 
tion  of  u>  for  a  beam  of  length  Z  -  50  ft  and  section  2  in.  x  4  in.  The 
material  characteristics  chosen  were  as  follows: 

p  -  3.0  density  relative  to  water 

tan  6  -0.1 

m 


Ef  -  27  x  106  psi 

Em  (ui)  “  Em  (o)  [1  +  i  lo8iou]*  1 0*  (sec-1) 

ER  (o)  -  1.35  x  10^  psi 
m 

v  -  0.4  v,  -  0.6. 

m  f 

The  variation  of  the  amplitude  with  u  is  shown  in  Figure  5  for 
elastic  vibrations  using  Efl(o)  as  elastic  modulus  and  for  viscoelastic 
vibrations.  Elastic  analysis  is  shown  by  a  full  curve,  and  viscoelastic 
analysis  is  shown  by  circles.  It  is  seen  that  for  the  first  elastic 
resonant  frequency  resonance  is  attained  for  the  viscoelastic  vibration 
also.  At  second  elastic  resonant  frequency,  there  is  already  a  small 
deviation  between  the  two,  but  still  the  viscoelastic  amplitude  is  very 
large . 
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Figure  5.  Amplitude  of  Lateral  Deflection  of  Beam, 


8.  SAMPLE  CALCULATIONS  OF  VISCOELASTIC  PROPERTIES 

Results  for  static  viscoelastic  properties  of  uniaxially  fiber 
reinforced  materials  have  been  given  in  Section  2.  It  has  been  seen 
that  considerable  simplification  is  achieved  when  it  is  assumed  that  the 
fibers  are  perfectly  rigid  relative  to  the  matrix.  However,  this  as¬ 
sumption  is  not  permissible  in  the  calculation  of  E^(t);  but  fortunately, 
the  assumption  is  not  needed  here,  since  (63)  is  an  explicit,  simple 
expression  incorporating  finite  elastic  stiffness  of  fibers. 

In  order  to  obtain  an  idea  of  the  approximation  involved  in 
the  assumption  of  rigid  fibers,  we  perform  a  sample  calculation  of  the 
axial  shear  relaxation  modulus  Ga(t)  for  the  case  of  elastic  fibers. 

The  Laplace  transform  of  Ga(t)  is  given  by  (59).  To  be  specific,  it  is 
assumed  that  the  matrix  behaves  according  to  th>?  simple  Maxwell  model. 
Thus,  the  deviation  matrix  stress-strain  relation  has  the  form 


(280) 


wheie  Gm  is  the  elastic  (initial)  shear  modulus,  nm  is  a  viscosity  coef¬ 
ficient,  and  a  dot  indicates  time  derivative.  The  theory  given  in  Ref¬ 
erence  [3]  Includes  the  use  of  time  differential  operator  stress-strain 
relations  of  phases,  such  as  (280),  for  computation  of  effective  visco¬ 
elastic  properties.  It  is  also  possible  "o  tic  in  (280)  directly  with  a 
more  general  stress-strain  relation  of  type  (49).  To  see  this,  take  the 
LT  of  (49) ,  which  is 

8lj(p)  -  2p  Gffl(p)  etj(p)  .  (281) 

Now  the  LT  of  (280)  is 


(282) 


(283) 


which  expresses  the  initial  elastic  response.  Taking  into  account  (283) 
and  (282)  in  (281)  and  rearranging  yields 


(P) 


2p% 

1  +  V 


Vp)  • 


(284) 
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where 


-  r  <285> 

m 

is  the  relaxation  time  of  the  Maxwell  matrix.  Comparison  of  (284)  with 
(281)  shows  that 


pCn(p) 


V 

1  +  V 


Now  (286)  Is  substituted  into  (59),  yielding 


(286) 


Ga  (p)  "  Gf  (1-c)  (1+T  p)  +  n  (l-c)p  '  1  +  T  p 
t  m  r.  m 


(287) 


Equation  (287)  Is  easily  Inverted  into  the  time  domain  by  the  method  of 
partial  fractions.  The  final  result  is 
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(288) 


where 


When  the  fibers  are  rigid,  $  -*•  «  and  then  (288)  reduces  to 


(289) 


G*(t) 
a 


Since  Gm  exp 


[-8 


^  i+c  fti 

■ - T^exp  "  • 

m  L  mJ 


(290) 


is  the  relaxation  modulus  of  the  Maxwell  matrix, 


(290)  is  in  accordance  with  the  general  result  (60). 


In  order  to  assess  the  numerical  importance  of  $,  (288)  has 
been  plotted  in  Figure  6  as  a  function  of  t/Tm  for  several  values  of 
<t>,  assuming  that  c  -  0.5.  It  is  seen  that  for  >  25 ,  which  is  a 
typical  value  for  fiber  reinforced  materials,  the  rigid  fiber  assumption 
gives  a  good  approximation.  Therefore,  this  Idealization  is  considered 
to  be  justified. 
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Matrix) 


Sample  calculations  of  viscoelastic  properties  of  fiber  rein¬ 
forced  materials  have  been  performed  for  a  polyisobutylene  matrix.  The 
matrix  properties  have  been  taken  from  Tobolsky  and  Castiff  [10].  Ini¬ 
tial  elastic  moduli  are 

Kffi(o)  -  2.32  x  105  psi,  Gffl(o)  -  1.85  x  105  psi. 

It  is  assumed  that  the  matrix  is  elastic  in  dilatation.  Plots  of  the 
shear  relaxation  modulus  G^t)  and  the  shear  creep  compliance  g^t)  are 
shown  in  Figures  7  and  8.  With  the  assumption  of  rigid  fibers,  the 
axial  shear  relaxation  modulus  G^(t)  is  given  by  (61).  Plots  of  G*(t) 
for  various  c  are  shown  in  Figure  9.  The  axial  shear  creep  compliance 
is  simply  given  by 

*>>  ■  £  '  (291> 

For  derivation  of  (291),  see  Reference  [3].  Equation  (291)  is  plotted 
for  various  c  in  Figure  10. 

The  plane  strain  bulk  relaxation  modulus  is  given  by  (58).  In 
the  present  case.,  it  is  assumed  that  the  matrix  is  elastic  in  dilata¬ 
tion. 

Therefore , 


Kffl(t)  -  Km  H(t)  ,  (292) 

where  H(t)  is  the  Heaviside  unit  step  function.  Then  (58)  assumes  the 
form 


K 

<(t>  ■  »<‘>  +  l(fc!y  G„«>  • 

A  plot  of  (293)  for  various  c  is  shown  in  figure  11. 

The  situation  with  respect  to  the  plane  strain  bulk  compliance 
k£(t)  is,  however,  not  so  simple.  The  relationship  between  Kj  and  k$  is 
expressed  in  terms  of  their  Laplace  transforms  by  (40) .  When  this 
expression  is  written  in  the  form 


kJ(p)  kj(p)  -  *  (294) 

P 

inversion  by  use  of  the  convolution  theorem  for  Laplace  transforms  leads 
to 
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Figure  7.  Shear 


a. 5 


igure  8.  Shear 


Figure  9.  Effective  Axial  Shear  Relaxation  Modulus 
G  ( t )  (Polyisobutylene  Matrix). 


Figure  10.  Effective  Axial  Shear  Creep  Compliance 
g*(t)  (Polyisobutylene  Matrix). 
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Figure  11.  Effective  Plane  Strain  Bulk  Relaxation  Modulus 
KT(t)  (Polyisobutylene  Matrix). 


t 


(295) 


K^(t-x)  k*(-r)  dt  -  t  . 

Equation  (295)  is  a  Volterra  integral  equation  with  kernel  K^(t-r)  and 
unknown  function  k$.  Introduction  of  (293)  into  (295)  yields 

K  ft  rt 

H(t-x)  k*(T )  dt  +  3 ( i-~c~)  Gnj(t-T )  k*(t)  dx  -  t  .  (296) 

o 

Now 

H(t-T)  ■  1  t> x  ,  (a) 

H(t-T)  -  o  t>T  .  (b)  (297'' 

It  is  seen  that  (297b)  is  impossible  since  o<r<t.  Therefore,  H(t-r)  in 
(296)  is  simply  replaced  by  unity.  Then  (296)  takes  the  form 

J‘  [h  +  itirsy  ]  k;(, )  d,  -  t  .  (298) 

The  integral  equation  (298)  has  been  solved  numerically  by  the 
method  of  Lee  and  Rogers  [11]  and  by  the  use  of  the  given  experimental 
information  K_  and  Gm(t) .  Plots  of  kJ  for  various  values  of  c  are  shown 
in  Figure  12. 

The  Young's  relaxation  modulus  E^(t)  is  given  by  the  simple 
expression  (63).  Here  it  is  not  permissible  to  assume  rigid  fibers, 
since  the  (63)  would  become  infinite.  While  (63)  is  easily  evaluated, 
there  is  little  point  in  doing  this;  for  it  has  been  shown  in  (63)  that 
for  the  usual  stiff  fibers  which  are  encountered  in  practice,  the  time 
dependence  of  (63)  is  negligible.  A  very  good  approximation  for  (63)  is 
simply 

*  % 

Eg  *=  cE^  -  constant,  t>o  .  (299) 

Similarly,  the  Young's  creep  compliance  e*(f)  -'.s  very  well  approximated 
by  the  inverse  of  (299).  Thus, 

*  %  1 

e  -  — =—  =  constant,  t>o  .  (300) 

B 

Another  important  viscoelastic  modulus  which  has  not  been  dis¬ 
cussed  in  detail,  so  far,  is  G^(t)  .  The  definition  of  this  modulus  is 
given  by  (23)  and  (25) .  Unfortunately,  the  composite  cylinder  assemblage 
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Figure  12.  Effective  Plane  Strain  Bulk  Creep  Compliance, 
kj(t)  (Polyisobutylene  Matrix). 


T 


1 


I 

I 


analysis  performed  in  Reference  [1]  has  not  yielded  a  closed  result  for 
the  corresponding  in  the  elastic  case.  Only  lower  and  upper  bounds 
have  been  obtained.  The  lower  bound  for  is  given  by 


m 


kE  +  2GE 
_ m _ m__ 

2(kE  +  GE) 
m  m 


(301) 


where  superscripts  E  denote  elastic  moduli;  see  Reference  [2],  page  19. 
The  upper  bound  for  Gj^  is  very  complicated  and  will  not  be  reproduced 
here.  In  the  event  that  the  bounds  are  close,  the  elastic-viscoelastic 
correspondence  principle  for  effective  properties,  as  expounded  in  Ref¬ 
erence  [3],  may  be  usgd  to  convert  (301)  into  the  approximate  Laplace 
transform,  G^(p)  of  Gj(t)  for  the  viscoelastic  fiber  reinforced  mate¬ 
rial.  (Strictly  speaking,  this  is  not  legal).  Considerable  simplifi¬ 
cation  is  achieved  if  the  fibers  are  taken  to  be  rigid  and  the  matrix 
is  taken  to  be  incompressible.  Then  it  is  easily  shown  that  (301)  in 
transform  domain  reduces  to  (60)  .  We  thus  obtain  the  approximation 


A  1 4-r-  A 

G  (t)  %  G  (t)  -  G  (t) 
T  i-c  m  a 


(302) 


It  is  to  be  expected  that  for  rigid  fibers  and  a  nearly  incompressible 
matrix  such  as  a  polymer,  (302)  will  provide  a  fair  approximation. 

We  now  turn  to  the  computation  of  effective  complex  moduli. 
Again,  a  polyisobutylene  matrix  is  chosen,  and  the  fibers  are  assumed  to 
be  rigid.  The  matrix  is  again  assumed  to  be  elastic  in  dilatation. 

This  implies  that  the  complex  bulk  modulus  (133)  of  the  matrix  satisfies 
the  conditions 


K  (oj)  -  K  *  constant,  (a) 
m  m 

K1  (a)  -  o  .  (b)  (303) 

m 

The  complex  shear  modulus  G„  (ico)  values  have  again  been  taken 
from  Reference  [10].  A  plot  of  Gjjj  (a;)  is  shown  in  Figure  15  (curve  c  * 
o),  and  a  plot  of  G^(ui)  is  g-'  fen  in  Figure  16  (curve  c  ■  o),  both  plots 
being  vs.  log  to. 

Kj1,  G^,  G**,  E*^,  and  E**  have  been  computed  on  the 
basis  of  Equations  (135),  (136),  (141),  (142),  (148),  and  (149), 
respectively.  The  results  for  various  c  are  shown  in  Figures  13  through 
18. 
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Figure  13.  Real  Part  of  Complex  Bulk  Modulus 
*R 

Kj  (oj)  (Polyisobutylene  Matrix). 
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Figure  14.  Imaginary  Part  of  Complex  Bulk  Modulus, 
Kj*(«)  (Polyisobutylene  Matrix). 
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Figure  15.  Real  Part  of  Complex  Axial  Shear  Modulus, 
G*R(w) (Polyisobutylene  Matrix). 
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Figure  17.  Real  Part  of  Complex  Axial  Young's  Modulus 
E,R(w)  (Polyisobutylene  Matrix). 


Figure  18.  Imaginary  Part  of  Axial  Young's  Modulus, 
E** (u>)  (Polyisobutylene  Matrix). 


9. 


LATERAL  VIBRATIONS  OF  FIBER  REINFORCED  VISCOELASTIC  PLATES 


It  has  been  shown  in  previous  sections  that  the  solution  of 
viscoelastic  problems  for  fiber  reinforced  materials  may  be  reduced  to 
the  solution  of  associated  elastic  problems  in  a  manner  identical  to 
viscoelastic  analysis  of  homogeneous  materials.  This  is  done  either  by 
application  of  the  elastic-viscoelastic  correspondence  principle  (Ref¬ 
erences  [12]  and  [13])  or,  for  steady-state  vibrations,  by  the  use  of 
suitably  defined  complex  moduli  (Reference  [14]).  In  this  section  we 
consider  the  lateral  vibrations  of  fiber  reinforced  viscoelastic  plates, 
for  which  the  associated  elastic  problems  are  those  of  lateral  vibra¬ 
tion  of  anisotropic  plates. 

DERIVATION  OF  GOVERNING  EQUATIONS 

The  governing  equations  for  anisotropic  plates,  based  on  the 
Kirchhoff  hypothesis,  are  easily  developed  by  integration  of  the  equa¬ 
tions  of  linear  elasticity  across  the  plate  thickness.  It  will  be 
remembered  that  the  displacements  are  of  the  form 

Ua  "  x3yc  <Vx2>  “  '  l-2  (a> 

u3  =  w  =  w  (x1,x2)  ,  (b)  (304) 

and  that  the  neglect  of  transverse  shear  deformation  leads  to  the 
relations 


Y 


a 


w, 


where 


w, 


3w 

3x 


(305) 


(306) 


The  equilibrium  equations  for  an  element  of  plate  are  easily  obtained, 
with  the  result 


3M 


aP 


?x. 


-  Q  =  o  . 


3Qa 

^  +  p  =  phw  , 

a 


(307) 


(308) 


where  dots  denote  partial  derivatives  with  respect  to  time.  It  should 
be  noted  that  the  rotatory  inertia  has  been  neglected  in  Equation  (307) 
This  is  consistent  with  the  neglect  of  transverse  shear  deformation  and 
will  limit  the  present  theory  to  frequencies  which  cannot  approach  the 


68 


frequency  of  the  first  thickness  shear  mode.  This  limitation  is,  how¬ 
ever,  not  severe  for  thin  plates. 


The  stress-strain  relation  for  anisotropic  plates  may  be 


written  as 


M 


'a  6 


BceByX  KyX  » 


(309) 


where 


K 


a  6 

* 


2  (Ya,B  +  t6,.) 


(310) 


The  plate  stiffness  tensor  B"  ,  is  defined  as 

*  aByX 

*  *  h3 

Ba8yX  "  BaByX  12  ’ 


(311) 


where  are  the  components  of  the  elasticity  tensor  for  plane  stress 

theory  for  the  material  referred  to  the  plate  axes.  It  should  be  noted 
that  could  also  be  interpreted  as  a  linear  integral  time  operator 

or  as  complex  moduli  for  either  formulation  of  the  viscoelastic  problem. 


Substitution  of  Equation  (310)  into  Equation  (309)  gives 

* 


M 


a  6 


-  B  .  ,  w,  .  . 

aByX  yX 


(312) 


The  governing  equation  for  anisotropic  plates  is  then  obtained 
from  Equations  (307),  (308),  and  (312);  the  result  is 


BaeV*  +  Cl,“  ■  1 


(313) 


where  p  is  the  transverse  load  per  unit  area. 

The  appropriate  boundary  conditions  are  that  on  each  edge  we 
must  prescribe  either  w  or  the  effective  shear  Va  and  either  the  slope 
or  the  bending  moment.  It  will  be  remembered  that  due  to  the  neglect  of 
transverse  shear  deformation,  the  problem  is  of  the  fourth  order  and 
does  not  allow  us  to  prescribe  the  shear  resultant  and  twisting  moment 
independently;  we  are  therefore  led  to  defining  the  effective  shear  Va 
as 


V 

a 


3M 


+  9x 


aB 


(no  sum  on  B) 


(314) 
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Expanding  Equation  (313),  we  obtain  the  governing  equation  in  the  form 


^4 

9  W 


Bllll  4  +  (B1211  +  B2111  +  B1112  +  B1121)  3 

3x^9x2 


+  (b2211  +  bj212 


4 

+  B1221  +  B1122  +  B2112  +  B2121)  2  2  +  (B2212  +  B2221  +  B1222 
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We  recall  that  the  stiffness  tensor  B*gy^  must  satisfy  the  symmetry 
relations  B*CyX  ■  B*aYA  “  Bagxy  “  B*yafi »  and  the  above  equation  there¬ 
fore  reduces  to 
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(316) 


Jf 

It  will  be  recalled  that  the  stiffness  components  Bag  .  appearing  in  all 
previous  equations  are  referred  to  plate  coordinates.  When  the  plate 
material  is  fiber  reinforced  with  unidirectionally  oriented  fibers,  the 
material  is  orthotropic  with  respect  to  axes,  which  are  respectively 
along  the  fiber  direction  and  perpendicular  to  the  fiber  direction  (Ref¬ 
erence  [1]).  It  is  clear  that  plate  stiffnesses  can  be  defined  with 
respect  to  material-oriented  axes  as 


(317) 


where  A^yA  are  the  components  of  the  elasticity  tensor  referred  to 
material  axes.  Clearly,  then,  the  stiffness  tensor  B0gy^  (or  Bagy^) 
obeys  the  transformation  law  for  a  surface  tensor  of  rank  4  which  is 
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where  aag  are  the  direction  cosines  of  the  plate  axes  with  respect  to 
the  material  axes.  If  we  consider  a  rectangular  plate  with  fibers 
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oriented  at  an  angle  6  to  the  direction,  it  is  easily  shown  that  the 
stiffness  tensor  referred  to  the  plate  axes  is  obtained  from  that  re¬ 
ferred  to  the  material  axes  by  the  matrix  equation 
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(319) 


Elements  of  the  matrix  [A]  are  given  as  follows: 
4 

A11  =  A22  *  COS  8 
4 

A12  =  A21  =  Sin  9 

2  2 

A13  =  A23  =  ^  s*n  0  COS  0 
2  2 

A..  =  A  .  =  2  sin  6  cos  6 
14  24 

2  2 

A31  =  A32  *  s*n  0  cos  0 

A  =  -  A 

33  13 

4  4 

A_.  =  sin  6  +  cos  9 

34 

A41  =  A42  =  A31 

A^  =  (sin^Q  -  cos^6)^ 


44 


14 


A,.^  =  -  A^  *  sin0  cos  0 

3 

A^2  =  ~  =  -  sin  0  cos0 

2  2 

A,.^  =  -  A^^  =  2  (sin  0  -  cos  0)  sin0  cos0 
2  2 

A,.^  =  -  A^  =  (sin  0  -  cos  0)  sin0  cos0 


(320) 
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It  is  clear  from  Equation  (319)  that  the  fiber  reinforced  plate  will,  in 
general,  be  anisotropic  for  0^0.  This  special  type  of  anisotropy  is 
often  denoted  by  the  term  related-orthotropic  [15]. 

VIBRATIONS  OF  ORTHOTROPIC  FIBER  REINFORCED  PLATES 


We  first  consider  the  simplified  case  of  fibers  oriented  along 
the  axis.  The  plate  is  then  orthotropic,  and  Equation  (316)  then 
reduces  to  the  form 
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where 
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(323) 


It  will  be  of  interest  to  consider  the  case  of  a  rectangular 
plate  of  length  "a"  and  width  "b"  simply  supported  along  all  edges  and 
subjected  to  a  sinusoidally  varying  concentrated  load  at  a  point 
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For  a  plate  simply  supported  along  all  edges,  the  boundary  conditions 
are 


w  =  =  0 

at  x^  «=  0, 

a  , 

(a) 

w  =  M22  *=  0 

at  x2  =  0, 

b  , 

(b) 

(324) 

and  we  note  that,  since  w  diminishes  along  these 
boundary  condition,  by  virtue  of  equations  (323) 

edges,  the  moment 
,  reduces  to 

S2vr  „ 

9x^ 

at  x^  =  0, 

a  , 

(a) 

a2 

~~§  ”  o  , 

'N _ ” 

at  x2  =  0, 

b  . 

(b) 

(325) 

If  we  consider  a  sinusoidally  time-varying  concentrated  load  of  magni¬ 
tude  P  applied  at  the  point  (£,,  °f  the  plate,  the  loading  function 
p(x^,  x^,  t)  may  be  expressed  as 


p(xx,  x2,  t)  =  P5  (Xj^-C-j^)  6  (x2~£2)  e'Ut  .  (326) 

It  Is  clear  from  the  boundary  conditions,  (324)  and  (325),  that  the  solu¬ 
tion  to  the  present  problem  may  be  obtained  as  a  Fourier  series  in  the 
mx^  nx2 

functions  sin  — g—  and  sin  ~ .  We  therefore  expand  p(x^,  x2,  t)  in 
such  a  Fourier  series  and  obtain 


,p  °°  00  mirf;  mTTX .  nnf,  nnx 

p(xl*  X2*  t)  =  ab  2  l  (sin  1“  Sin  ~  Sln  ~  Sin  ~) 

m=l  n=l 

We  now  seek  a  solution  for  w(x^,  x2,  t)  in  the  form 


e 


lOJt 


(327) 


00  ®  mTTX  nirx_ 

w(x  ,  X  ,  t)  *  [  l  q  (t)  sin  — —  sin  — —  .  (328) 

~  i  *  mn  3  D 

m=i  n=l 

Substituting  Equations  (327)  and  (328)  into  Equation  (321)  and  equaling 
coefficients  of  like  terms  in  the  Fourier  series  yields  one  ordinary 
differential  equation  for  the  functions  qmn(t);  the  result  is 
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(329) 
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This  equation  may  be  written  in  the  form 
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where 
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It  will  be  of  particular  interest  to  consider  only  the  steady-state 
solution  of  Equation  (330),  since  in  the  viscoelastic  case  transients 
will  be  damped  out.  Thus,  we  consider  a  solution  of  the  form 


V,(t> 


\U)t 


mn 


and  substitute  this  expression  into  Equation  (330).  We  then  obtain  the 
steady-state  solution  as 


^n(t)  " 

mn 
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nut 
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(333) 


Thus,  the  elastic  solution  to  the  problem  may  be  written  in  the  form 


lut 

00  00  p  e  m7Tx1  n7Tx„ 

w(x1,  x?,  t)  -  l  l  - 2  sin  ~~a —  sln  ~b —  *  (334) 

m-1  n-1  u)  -  uj 
mn 

As  was  clear  in  the  case  of  the  beam  solution  of  Section  7,  this  solu¬ 
tion  may  be  written  in  the  form 

00  00 

w(xlf  x2,  t)  -  elwt  l  l  W^  (xx,  x2)  ,  (335) 

m-1  n-1 
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where 
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(336) 


The  passage  to  the  viscoelastic  solution  may  now  be  accomplished  in  a 
manner  entirely  analogous  to  that  of  Section  7.  Thus,  for  the  visco¬ 
elastic  case,  we  introduce  into  Equation  (331)  appropriate  complex  stiff¬ 
nesses  for  D*m,  D*222»  anc*  ^1122*  These  complex  stiffnesses  are  com¬ 
pletely  defined  by  Equation  (331)  if  the  components  or  the  elasticity 
tensor  are  replaced  by  the  corresponding  complex  moduli.  They  can  be 
written  explicitly  in  terms  of  complex  moduli  as 
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expression 


Under  these  conditions,  Equation  (331)  is  then  replaced  by  the 
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(338) 
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It  should  be  noted  that  wmn  and  6mn  are  both  functions  of  frequency  u>, 
since  the  complex  moduli  and  stiffnesses  are  frequency  dependent.  Thus, 
the  viscoelastic  solutions  to  the  problem  can  now  be  written  in  the  form 

QO  CXI 

w^,  x2,  t)  «  elut  l  l  W™  (x.^  x2)  ,  (339) 

m*l  n-1 
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Equation  (340)  may  be  written  as 
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Following  the  development  previously  presented  in  Equations  (260)  through 
(263),  we  can  obtain  the  dynamic  viscoelastic  solution  to  the  problem. 

To  illustrate  the  dynamic  viscoelastic  solution  (339),  we  con¬ 
sider  a  rectangular  plate  of  length  "a"  and  width  "b"  simply  supported 
along  all  edges  and  subjected  to  a  sinusoidally  varying  concentrated 

load  at  the  center  fj,  .  See  Figure  19.  The  dimensions  of  the  plate 
are 

a  -  60  in.,  b  =  30  in.,  h  ■  0.375  in. 

The  matrix  material  of  the  plate  is  taken  to  be  polyisobutylene,  whose 
properties  are  given  in  Section  8.  It  is  assumed  that  the  matrix  is 
elastic  in  dilatation.  The  material  properties  of  the  fiber  and  the 
density  of  the  plate  are  as  follows: 

=  5.83  x  10^  psi 


Gf  =  4,375  x  10 


psi 
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p  =  3.0  (density  relative  to  water) 
c  =  0.5  (volume  fraction  of  fiber). 


From  (339)  and  (340),  the  deflection  of  the  plate  at 
is  found  to  be 
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where 
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Using  (337),  the  deflection  at  the  center  of  the  plate  can 
be  computed  from  (342).  The  variations  of  the  amplitude  with  u>  are 
plotted  in  Figure  20.  As  in  the  case  of  beams,  it  is  seen  that  there  is 
no  true  resonance.  However,  for  forcing  frequencies  close  to  the  first 
natural  frequency  of  the  plate,  large  response  amplitude  is  encountered 
and  resonance  is  attained. 

Finally,  we  note  that  the  algebraic  and  numerical  calculations 
for  this  case  are  at  least  one  order  of  magnitude  more  complex  than 
those  for  the  case  of  beams,  which  was  treated  previously.  This  is  due 
to  the  appearance  of  four  distinct  moduli  in  Equation  (331),  each  of 
which  has  different  frequency  dependent  real  and  imaginary  parts;  this 
causes  the  "loss  tangent",  tan  6mn,  to  depend  on  the  mode  under  con¬ 
sideration.  Some  simplifications  may,  however,  be  possible  for 
particular  materials. 
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AMPLITUDE 


10.  ELEMENTARY  THEORY  OF  BENDING  AND  SHEARING 

or  viscoelastic  rim  rein'forcfd_beam$ 


We  consider  a  beam  of  constant  cross  section.  The  system  of 
axes  to  which  the  beam  is  referred  is  shown  in  Figure  21.  It  is  assumed 
for  the  sake  of  simplicity  that  the  y-axis  is  an  axis  of  symmetry  of  the 
cross  section.  The  beam  is  uniaxially  reinforced  by  fibers  in  the  x- 
direction. 


It  will  be  our  purpose  to  develop  a  theory  of  such  beams  which 
is  analogous  to  the  Euler-Bernoulli  theory  of  homogeneous  beams.  We 
note  in  this  respect  that  for  homogeneous  elastic  and  viscoelastic  beams 
in  pure  bending,  the  Euler-Bernoulli  theory  is  rigorous,  provided  that 
the  bending  moments  on  the  terminal  beam  sections  are  .pplied  through  a 
normal  stress  distribution  which  varies  linearly  over  the  section.  Sup¬ 
pose  that  such  a  stress  distribution  is  applied  on  the  terminal  sections 
of  a  uniaxially  fiber  reinforced  beam.  The  beam  may  be  thought  of  as 
being  represented  by  the  composite  cylinder  assemblage  model.  Since 
there  are  a  very  large  number  of  composite  cylinders  passing  through  the 
section,  the  normal  stress  oxx  on  any  composite  cylinder  end  section  is 
uniform.  Then  it  is  known  that  in  such  a  composite  cylinder,  the  axial 
strain  cxx  will  be  uniform  throughout  fiber  and  matrix  (see  Reference 
[1]),  except  for  some  local  perturbation  at  the  end  regions.  Since  the 
strain  cxx  is  proportional  to  the  stress  csxx,  it  follows  that  cxx  is 
very  nearly  linearly  distributed  over  the  section  of  the  cylinder.  We 
therefore  conclude  that  in  pure  bending,  the  Euler-Bernoulli  theory 
should  be  applicable  with  high  accuracy  for  a  fiber  reinforced  beam. 

Because  of  the  assumption  that  plane  sections  remain  plane  and 
since  y  is  an  axis  of  symmetry  of  the  section,  we  write,  as  in  the  usual 
strength  of  materials  development, 


Sex  (x'!,>t) 


R(x,t)  ’ 


(344) 


where  R(x,t)  is  the  radius  of  curvature  of  the  deformed  beam,  which  here 
depends  also  on  time  because  of  the  time  dependence  of  the  material. 

The  average  stress  oxx  on  the  section  of  a  composite  cylinder  is  then 
given  by 


5xlI(x,y.O 


9  Ev*(X’y,T) 

(t-r)  — ^ -  dr 


(345) 


Note  that  the  y-coordinate  represents  positions  of  central  axes  of  com¬ 
posite  cylinders  only.  Since  there  are  very  many  of  these,  y  may  be 
approximated  by  the  continuous  variable  y. 

Introducing  (344)  into  (345),  we  have 
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In  the  absence  of  axial  forces  on  the  beam,  equilibrium  demands  that 


°vv(x*y»t)  dA  -  0  .  (348) 

J  (A)  X 

Introducing  (346)  into  (348)  ,  it  follows  that 

ydA  =  0  .  (349) 

J(A) 

This  implies  that  the  linear  distribution  (346)  must  pass  through  the 
centroid  of  the  section  and  that  the  stress  and  strain  vanish  at  the 
centroid . 

From  moment  equilibrium  at  any  section,  we  have 

M  (x,t)  =  M(x,t)  =  o  (x,y,t)  ydA  .  (350) 

j  (A)  X 

Introducing  (346)  into  (350),  we  find  that 
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M(x,t) 


B(x,t)  I  , 


(351) 


where 


I  -  I 


zz 


y2dA 


(352) 


Is  the  moment  of  Inertia  with  respect  to  the  z-axis.  Combining  (351) 
with  (346) ,  we  have 

°xx(x,y,t)  "  ^  y  ’  (353) 

which  is  the  exact  analogue  of  the  strength  of  materials  stress  for  an 
elastic  beam.  It  is  thus  seen  that  the  stress  oxx  is  independent  of  the 
viscoelastic  nature  of  the  material.  It  should  be  noted  that  in  accor¬ 
dance  with  the  present  development,  the  moment  M  is  reckoned  to  be 
positive  if  it  produces  tension  below  the  neutral  axis  and  compression 
above  the  neutral  axis  of  the  section. 

In  order  to  analyze  the  deflection  of  the  center  line  w(x,t), 
we  assume  this  deflection  to  be  small  in  comparison  to  section  dimen¬ 
sion.  Then 


R(x,t)  £  -  3  .  (354) 

3x 

The  negative  sign  in  (354)  stems  from  the  choice  of  positive  y  downward 
and  bending  moment  sign  convention.  It  follows  from  (354)  and  (347) 
that 
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dx  . 


Introducing  (355)  into  (351) ,  we  have 
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(356) 


Equation  (356)  is  the  governing  equation  for  the  beam  deflection.  End 
conditiona  for  w  have,  of  course,  to  be  specified. 

The  present  form  of  (356)  is  somewhat  inconvenient.  We  shall 
now  proceed  to  transform  it  into  a  different  form.  For  this  purpose  let 
us  take  the  time  Laplace  transform  of  both  sides.  Then 
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M(x,p)  «  -  I  — ^  p  Efl  (p)  w  (x,p)  ,  (357) 

where  we  have  used  the  notation  (3)  for  Laplace  transforms.  Let  us  con¬ 
sider  the  very  important  case  where  the  load  is  applied  to  the  beam  at 
time  t  -  0  and  is  then  held  constant. 

In  that  case, 


M(x,t)  -  M(x)  H(t)  ,  (358) 

where  H(t)  is  the  Heaviside  unit  step  function  defined  as 

rO  t  <0 
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It  follows  that 


From  (45), 


M(x,p)  -  -  M(x)  . 


(359) 
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pefl(p) 


where  efl(p)  is  the  transform  of  the  Young's  creep  compliance  efl(t)  which 
enters  into  the  relation 


e  (t) 
xx 


.  3o  (t) 

e  (t-x)  — r -  dx  . 

a  ox 


(360) 


Thus,  e  (t)  is  the  longitudinal  strain  produced  in  a  fiber  reinforced 


- ,  -o'-'  —  -  — O -  - r 

material  by  constant  unit  average  stress 


;xx* 


Introducing  (359)  and  (45)  into  (357),  we  find  that 


3  w(x.d)  M(x)  -*  ,  . 

3V  ■  -  i  %  <p>  • 


(361) 


which  can  be  directly  inverted  into  the  time  domain  to  give 


2-aifSl  .  .  i!M.  e*  (t) 

3x  1  3 


(362) 
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Equation  (362)  is  the  analogue  of  the  well-known  differential  equation 
for  an  elastic  beam.  As  is  usually  done  for  elastic  beams  we  shall  now 
make  the  assumption  that,  although  derived  for  the  case  of  pure  bending, 
Equation  (362)  is  also  applicable  when  shear  forces  are  present.  Then 
(362)  represents  the  differential  equation  of  viscoelastic  fiber  rein¬ 
forced  beams  in  bending. 

We  now  consider  end  conditions  for  w.  Typical  end  conditions 

are 


simple  support  at  xq 

w(xQ,t)  -  0  (a) 

built-in  support  at  x^ 

&  <V£)  ■  0  •  (b) 

(363) 

Set 

w(x,t)  -  W(x) 

Then  Equation  (362)  reduces  to 

* 

e  (t)  . 
a 

(364) 

.  _  Mix! 

.2  I 

dx 

(365) 

Also,  introduction  of  (364)  into  (363) 

gives 

W(x  )  -  0 
o 

(a) 

dW(x  ) 

-T-°_  _  o  . 
dx 

(b) 

(366) 

It  is  seen  that  (365)  together  with  (366)  represents  an  elas¬ 
tic  deflection  problem  with  unit  Young's  modulus.  The  viscoelastic 
solution  is  then  obtained  from  (364).  Since  in  any  elastic  beam  deflec¬ 
tion  problem  El  occurs  in  the  deflection  denominator,  we  can  now  formu¬ 
late  the  following  simple  rule:  In  order  to  find  the  deflection  of  a 

viscoelastic  beam ,  simply  take  the  elastic  deflection  and  replace  the 
reciprocal  of  the  Young's  modulus  by  the  viscoelastic  Young '8  compli¬ 
ance. 


As  an  application,  consider  a  simply  supported  beam,  of  length 
l,  which  is  loaded  uniformly  by  p0  per  unit  length.  The  elastic  deflec¬ 
tion  is 


wE(x)  -  (Jt3-2£x2+x3) 

24E  I 
a 


(367) 
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The  viscoelastic  deflection  is 


t.r.  P  (O  q  9  O 

wVE(x,t)  =  -  U3-2£x2+x3)  .  (368) 

The  maximum  deflection  is  at  the  center  x  =  1/2 


wVE(£/2,t)  =  e*(t)  .  (369) 

It  has  been  shown  in  Reference  [3]  that  the  viscoelastic  effect 
in  uniaxial  stress  is  negligible  for  conventional  fiber  reinforced  mate¬ 
rials  in  which  the  fibers  are  very  stiff  in  comparison  to  the  matrix. 

In  other  words,  if  the  creep  compliance  ea(t)  is  replaced  by  the  Initial 
compliance  e*(o),  the  error  produced  is  insignificant;  but  according  to 
the  development  given  in  Reference  [3], 

e*(°)  =  ,  (370) 

E  (o) 
a 

where  from  (63) 

E  (o)  =  (1-c)  E  (o)  +  cE  ,  t^.o  .  (371) 

a  mi 

(It  should  be  recalled  that  c  is  the  volume  fraction  of  fibers  and  that 
Em(o)  is  the  initial  (elastic)  Young's  modulus  of  the  matrix.)  Thus, 
(368)  assumes  the  form 

wVE(x,t)  =  — - -  (£3-2£x3+x3)  .  (372) 

24E  (o)I 
a 

Equation  (372)  is  no  logger  time  dependent  and  is  precisely 
the  elastic  solution  (367)  with  Ea(o)  used  as  an  elastic  Young's  mod¬ 
ulus.  It  is  also  quite  clear  that  similar  conclusions  would  hold  for 
bending  deflection  analysis  for  any  loading,  and  we  may  thus  state  the 
following  important  conclusion:  Under  time-constant  load,  the  bending 
deflections  of  uniaxially  fiber  reinforced  viscoelastic  beams,  in  which 
the  fibers  are  very  much  stiffer  than  the  matrix,  may  be  approximated 
with  high  accuracy  by  the  time-independent  deflection  of  elastic  beams 
with  Young's  modulus  taken  as  the  initial  E*(o). 

We  now  turn  to  the  effect  of  shearing  forces  on  the  deflection 
of  fiber  reinforced  beams,  and  we  shall  adopt  the  usual  strength  of 
materials  approach.  It  will  first  be  recalled  th3t  bending  stresses  are 
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given  by  (353).  In  the  event  of  constant  load  for  t£o,  the  moment  is 
given  by  (358);  consequently,  the  bending  stresses  are  given  by 


5xx(x,y)  =  y’  t'°  ‘  (373) 

It  will  be  recalled  that  the  shear  stress  oXy  follows  accord¬ 
ing  to  strength  of  materials  simply  from  equilibrium  considerations  with 
the  use  of  the  bending  stresses  only.  Since  (373)  is  the  exact  analogue 
of  the  strength  of  materials  bending  stress,  it  follows  immediately  that 
the  shear  stress  oXy  is  also  given  by  the  strength  of  materials  expres¬ 
sion.  Thus, 


5xy(lt,y> 


S  (x)Q(y) 
b  (y)  I 


(374) 


where 


S(x)  =  Shear  force  at  x-section, 

Q(y)  =  First  area  moment  with  respect  to  the  neutral 
axis  of  the  part  of  the  section  above  y, 

b (y )  =  Width  of  the  section  at  height  y, 

I  =  Moment  of  inertia  I 

zz 

The  overbar  in  (374)  denotes  local  averages.  The  maximum  shear 
stress  always  occurs  at  y  =  0.  It  may  be  conveniently  written  in  the 
form 


c  „  =  °  (x>°)  =  k  » 


xy  max 


xy 


SfcLl 

A 


(375) 


where  A  is  the  area  of  the  section  and  k  is  a  geometrical  section  shape 
factor  given  by 


k  _  (376) 

lb (o )  •  K  7 

3  4 

For  a  rectangular  section,  k  =  — ;  for  a  circular  section,  k  =  —  . 

The  shear  strain  produced  on  the  neutral  axis  by  the  stress 
(375)  is  given  by 


Scy(x,o)  =  2  8I(t)  ^xy(x>0)  ’  (377) 

This  expression  needs  some  explanation.  Since  the  stress  (375)  is  con¬ 
stant  in  time  for  t;o,  the  strain  produced  is  directly  proportional  to 
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the  appropriate  creep  compliance.  This  creep  compliance  is  g*(t)  since 
the  shearing  action  proceeds  in  places  parallel  and  transverse  to  the 
fibers . 


On  the  basis  of  the  treatment  given  in  Section  9,  the  creep 
compliance  for  stiff  fibers  can  be  very  well  approximated  by  the  creep 
compliance  for  rigid  fibers.  Accordingly,  we  use 

<(t)  *  £  g„<‘>  '  <378> 

According  to  strength  of  materials  theory,  there  is  an  addi¬ 
tional  shear  deflection  ws  which  is  superposed  on  the  usual  Euler- 
Bernoulli  bending  deflection  w^.  Since  w^  produces  no  shear,  we  have 

1  3w 

"  (x,o)  =  4  T-S-  .  (379) 

xy  2  3x 

Substituting  (375)  into  (377)  and  substituting  the  resulting  expression 
into  (379),  we  find  that 

9w 

if  -  A  S(x)  ga(t)  •  (380) 

Now  from  equilibrium 

=  -  p(x)  ,  (381) 

where  p(x)  is  the  load  per  unit  length  of  the  beam,  taken  as  positive  in 
the  y  (downward)  direction.  Differentiating  both  sides  of  (380)  and 
using  (381),  we  obtain 


,2 

<3  w 


3x 


f  =-  Jp(x)ga(t) 


(382) 


This  is  the  governing  equation  for  the  shear  deflection  of  the  beam. 

Equation  (382)  is  similar  to  (362)  and  can  be  solved  by  analo¬ 
gous  methods.  We  set 


w  (x , t )  =  W  (x)  g  (t)  . 

S  S3 


(383) 


Introducing  (363)  into  (382),  we  obtain  for  Wg  the  differential  equation 


.2 

d  w 


dx 


s  k  .  . 

- A 


(38A) 
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As  an  example,  we  consider  again  the  case  of  a  simply  sup¬ 
ported  uniformly  loaded  beam  with  p0  per  unit  length  for  t^o.  The  bend¬ 
ing  moment  is  then  given  by 


P  x 

M(x)  =  — ~  (£-x)  . 


(385) 


Introducing  (385)  into  (365)  and  integrating,  we  have 


Vx) 


_  Zo  (ixi  _  +  cb  +  cb 

21  l  6  12  )  L1  +  z 


(386) 


From  (364)  and  (386) , 


,  ,  poea(t)  /  tx3  x4\  rb  *,  . 

Vx't)  -  -  ~2I -  -  12  +  Clea(t)  x 


b  * 

C?e  (t) 
Z  a 


(387) 


Integration  of  (384)  for  constant  p^  yields 


Hs(x)  -  -  2r*2  +  cix  +  c2  • 


(388) 


Consequently,  from  (383), 


ws(x,t)  =  - 


kpoBa(t)  2 


2A 


x  +  C  g  (t)x  +  C  g  (t) 
la  z  a 


(389) 


The  total  deflection  is  then  given  by 


w(x,t)  =  w  (x,t)  +  w  (x,t) 
D  S 


p  e  (t)  /  3 

*o  a  '  °'r 


21 


f«x i  S\  ^ 

\  6  12 


«l(t)  2 

2A - X  +  Ci(t)  x  +  c2(t)  ,  (390) 


where 


Cl(t)  =  C1  ea(t)  ■'  C1  8I(t)  *  (a) 


C,(t)  =  C*  e*  (t)  +  C®  g*(t)  .  (b) 

2  2  a  2  a 


(391) 


The  end  conditions  of  simple  support  are 
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w(o ,  t)  =  w(£ , t)  =  o  . 


(392) 


Insertion  of  (390)  into  (392)  yields 


C2(t)  =  o  , 


(393) 


p  £  *  p  k£  ^ 

Cl(t)  "  241  ea(t)  +  2A  8a(t)  * 


(394) 


Hence,  finally, 


p  x  o  o  *  P  kx 

w(x,t)  =  241  (Z  -  Ux  +  x  >  ea(t)  +  ~2A" 


(£-x)  g  ft)  .(395) 

a 


The  first  part  of  (395)  is  recognized  as  the  bending  deflec¬ 
tion  wfo  (368),  while  the  second  part  is  the  shear  deflection  wg. 

Inspection  of  the  expressions  obtained  for  w^  and  wg  shows  a 
fundamental  difference.  According  to  previous  discussion  the  time  de¬ 
pendence  of  wjj  is  negligible,  and  for  all  practical  purposes  this  is  an 
elastic  deflection  which  does  not  vary  in  time.  However,  g*(t)  in  wg 
has  a  considerable  time  variation.  Indeed,  the  rigid  fiber  approximation 
(378)  shows  that  ga(t)  is  directly  proportional  to  the  shear  compliance 
of  the  matrix  gm(t) ,  and  this  compliance  can  assume  very  large  values, 
and,  theoretically  speaking,  can  even  become  unbounded.  We  thus  arrive 
at  the  conclusion  that  a  viscoelastic  fiber  reinforced  beam  can  become 
unserviceable  due  to  excessive  shear  deflection. 

As  an  example,  consider  the  deflection  (395)  at  the  center  of 
the  beam,  where  it  attains  its  maximum.  We  have 


w(H/2 , t) 


5p./  ,  kV2  l-c  ,M 

384E*I  8A  1+C  m 
a 


(396) 


where  we  have  neglected  the  time  dependence  of  eg(t)  in  accordance  with 
previous  discussion  and  have  replaced  it  by  1/E*.  Also,  the  rigid  fiber 
approximation  (378)  has  been  used  for  ga(t).  Denoting  the  first  term 
on  the  right  side  of  (396)  by  6e,  we  rearrange  (396)  in  the  form 


w(£/2 , t) 


6 

e 


2  * 
48kp  E 
_ a 


5£ 


l-c 

1+c 


where  .■>  is  the  radius  of  gyration  of  the  section;  i.e., 


(397) 
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p  -J"y.  (398) 

A  plot  of  (397)  is  shown  in  Figure  22  for  a  beam  of  rectangu¬ 
lar  section.  The  dimensions  of  the  beam  are 

b  =  2.0  in. 

h  =  4.0  in. 

I  =  60.0  in. 

The  material  characteristics  chosen  are  as  follows: 

=  27.0  x  106  psi 

E  (o)  =  3.62  x  105  psi 
m  r 


The  matrix  is  taken  to  be  polymethyl  methacrylate  (Reference  [9]), 
whose  shear  creep  compliance  is  shown  in  Figure  23. 


The  deflection  is  plotted  against  logarithmic  time  starting 
out  at  1  sec,  when  the  shear  deflection  is  negligible  in  comparison  to 
the  bending  deflection.  The  initial  deflection  doubles  (assuming  con¬ 
tinued  linearity  in  log  t)  after  10^7  hrs.  The  enormous  amount  of  time 
required  to  double  the  deflection  shows  that  there  is  no  danger  of 
failure  within  a  reasonable  lifetime  of  the  beam. 
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6  £1 


Figure  22.  Lateral  Deflection  of  Beam 


08  0 


hear  Creep  Compliance  of  Polymethyl  Methacrylate. 


11.  TORSION  OF  FIBER  REINFORCED  CYLINDER 


The  theory  of  torsion  of  uniaxially  reinforced  cylinders  when 
the  fibers  are  in  the  direction  of  the  generators  fortunately  proves  to 
be  very  simple.  Consider  a  cylinder  whose  x^-axis  is  in  the  generator 
and  fiber  direction,  while  X2,  X3  are  the  transverse  axes  in  the  plane 
of  the  transverse  section  A.  We  make  an  isotropic  elastic  torsion-type 
assumption 


=  ct(t)  4>  (x2,  x3)  , 

(a) 

=  -  a(t)  x3,  x1  , 

(b) 

=  u(t)  x2,  x3  , 

(c) 

(399) 

where  u(t)  is  the  time-variable  angle  of  twist  per  unit  length  of  cylin¬ 
der  and  <t>  is  the  torsion  function.  The  only  nonvanishing  strains  are 
then 


€12  =  2 


i£_  _  x 

9x2  3 


E13  “  2  a(t)  3x^  +  X2 


,  (a) 


.  (b) 


(400) 


Because  of  (20),  (21),  and  (24),  the  only  nonvanishing  stresses  are  then 

ft 
J  o 
ft 


°12  =  2 


*  (T) 

Ga(t-T)  ~T7~dT  *  (a) 


°13  ‘  2 


*  3e13<t) 
Ga(t-T>  3t  dT  ‘ 


(401) 


Inserting  (400)  into  (401)  we  find 

rt 


0.,  =  _ 

12  (  dx2  3 


G*(t-i)  daiLl  ^  ^  (a) 


13 


|^-+  x2  j  j  G*(t-t)  di  .  (b) 

'  1  J  rt 


(402) 


93 


Noting  that 


L(t) 


u  * 

G  (t-T) 
a 
o 


da  (t  ) 
dx 


dx 


(403) 


where  L(t)  is  a  known  function  of  time,  the  only  surviving  equilibrium 
equation  is 


3o 


12 


3o 


13 


3x, 


3x, 


=  o . 


(404) 


*2  va^ 

Substitution  of  (402)  into  (404)  yields  the  equation 


V2<(> 


2  2 

3  $  +  3  if 


3x 


2 


o  . 


(405) 


3x, 


Satisfaction  of  the  zero  traction  boundary  conditions  on  the 
curved  surface  of  the  cylinder,  for  the  whole  time  range,  leads 
precisely  as  in  the  elastic  torsion  problem  (see,  e.g.,  Sokolnikoff  [7]) 
to  the  <}>  boundary  condition 

^  *  x3  cos  (x2,  n)  -  x2  cos  (x3,  n)  ,  (406) 

on  the  contour  C  of  the  section  A,  where  n  is  the  normal  to  the  contour. 
Equations  (405)  and  (406)  define  4>  uniquely  except  for  an  immaterial 
arbitrary  constant.  We  conclude  that  the  present  function  $  is  the  same 
as  that  in  the  isotropic  elastic  torsion  problm. 

The  torque  on  the  cylinder  is  given  by 


'X2°13 


X3°12 


)  dx„  dx. 


(407) 


Substituting 


Ml- 

We  note  that 
given  by 


(402)  with  (403)  into  (407),  we  find  that 


L(t) 


[ 

, 

x22  + 


’■32  +  *2  Ifc  -  *3  Hr)  dx2  dx3  •  (408) 


the  initial  elastic  torsional  rigidity  of  the  cylinder  is 


el. 


=  G* (o) 
a 


+  x„  +  x 


_ _  1L_ 


2  3x, 


3  3x, 


dx2  dx^. (409) 
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Therefore  (408)  may  be  rewritten  as 


Hj(t)  -  *4^-  D_„  . 


it  0 

G  (o)  6U 
a 


(410) 


Furthermore,  the  stresses  (402)  in  the  cylinder  can  be  written  as 


_kllL 


'12  * 

1  a(o)G  (o) 
a 


-  x3  |  o(o)Ga(o)  ,  (a) 


Jd&L 


13  * 

a (o)G  (o) 
a 


ii-  +  x 

3x3  2 


a(o)G  (o)  .  (b) 

a 


(411) 


It  is  seen  that  if  "elastic"  stresses  are  defined  by 


efc. 


12 


M _ x 

3x2  3 


a(o)Ga(o)  ,  (a) 


°12 


el. 


ii_  +  x 

9x3  2 


a(o)Ga(o)  .  (b) 


(412) 


then  (411)  assumes  the  form 


el. 


12  *  12 

°(o)G  (o) 

a 


,  (a) 


L(t) 


el. 


13  *  "13 

a(o)G  (o) 
a 


.  (b) 


(413) 


As  an  application  let  us  consider  the  case  when  a  viscoelastic 
fiber  reinforced  cylinder  is  subjected  to  a  constant  torque  Mp0  for 
t  >  o.  In  that  event 


Hj(t)  -  Hp0  H(t) 


(414) 


where  H(t)  is  the  Heaviside  unit  step  function.  How  take  the  Laplace 
transforms  of  (403),  (410),  and  (414),  and  combine  the  results.  This 
gives 
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(415) 


V  Ga(o) 

*(p>  =  -§^ 


P  Ga(P)De£. 


If  (29)  is  recalled,  it  is  seen  that  (415)  assumes  the  form 

* 


a(p)  = 


<=a(o) 


el. 


ga(p)  • 


(416) 


This  result  can  be  directly  inverted  into  the  time  domain,  obtaining 

* 


a(t) 


V  Ga(o)  * 


D 


el. 


g,(t)  . 


(417) 


To  find  the  stresses,  take  the  Laplace  transform  of  (413), 
using  (403),  The  result  is 


P  Ga(p)a(p) 

- -  a  ,  (a) 

a(o)G  (o)  ^ 

a 


A  £  A 

P  G  (p)a(p) 


13 


a(o)G  (o) 
a 


el.  „  \ 

*—  °13  ‘  (b) 


(418) 


Introduce  (415)  into  (418)  to  obtain 


el.  V 

0„„  =  _■  /  v  V -  ,  (a) 


12  12  p  a (o)  D 


el. 


el.  V 

a,„  *  a,„  — r~x— r -  .  (b) 


13  13  p  a (o)  D 


el. 


(419) 


Because  of  (417), 


V  Ga(0)  * 


a(0)  =  —  —  - -  g_(o)  , 

el.  3 


(420) 
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but 
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Furthermore,  it  has  been  shown  in  Reference  [3]  that 


“  TT§  en(t)  *  (427> 

where  gra(t)  is  the  shear  creep  compliance  of  the  matrix.  Inserting 
(426)  and  (427)  into  (417),  we  find  that 


a(t) 


V  i- 


4  1+c  8m(t) 


(428) 


It  is  seen  that  if  gm(t)  becomes  unbounded,  which  is  quite  often  the 
case,  then  a(t)  becomes  unbounded. 
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12.  QUASI-STATIC  DEFORMATION  OF  VISCOELASTIC  FIBER  REINFORCED  PLATfcS 


We  consider  here  the  case  of  quasi-static  bending  by  trans¬ 
verse  loads  of  a  thin,  rectangular  fiber  reinforced  plate.  The  rein¬ 
forcement  direction  is  taken  as  x^.  The  middle  surface  of  the  plate  of 
thickness  h  is  referred  to  the  plane  Cartesian  system  of  axes  x^,  X2. 
The  X3~axis  is  chosen  normal  to  the  middle  plane,  pointing  dovmwaid 
(see  Figure  19) . 

The  deflection  of  the  middle  surface  is 

u^(x^,  x^,  o,  t)  =  w(x^,  x2 ,  t)  ,  (429 ) 

where  t  is  the  time.  We  use  the  usual  Kirchhoff  hypothesis  of  thin 
plate  bending,  according  to  which  the  deformations  are  given  by 

3  w 

U1<V  x2’  V  c)  -  -  *3  IT  ■  (',30) 


U2(V  V  V  ')  '  -  x3  3^  '  (431) 

The  strains  in  x^  -  X2  planes  associated  with  (429)  ,  (430), 
and  (431)  are  then  given  by 


e,  =  -  x 


3  w 


11 


3  .  2  ’ 

3X1 


E22  =  '  X 


,2 

3  w 


3  .  2  ’ 

3x2 


E12  "  X 


a2 

3  w 


3  3x^3x2 


(432) 


(433) 


(434) 


It  will  be  convenient  to  continue  the  analysis  in  Laplace 
transform  domain.  The  Laplace  transforms  of  (432)  through  (434)  are 
written  compactly  as 


ij  ^xl*  X2 ’  x3 ’  “  -  X' 


32w(xlf  x2,  p) 
3x^  9Xj 


(435) 


i.J  =  1.2 
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The  Laplace  transformed  viscoelastic  stress-strain  law  of  the 
uniaxially  reinforced  plate  material  may  be  written  in.  the  form 


pE  (p) 

d 


Op) 

pEa(p) 


%(P>  '  x  °22 

+  -*  » 
pEa(p)  pET(p) 


(436) 


(437) 


12 


12 


2pG  (p) 


(438) 


Here  Ea(p)  is  the  LT  of  the  axial  Young's  relaxation  modulus, 

££(p)  is  the  LT  of  the  transverse  (normal  to  fibers)  Young's  relaxation 
modulus,  G*(p)  is  the  LT  oi  the  axial  shear  relaxation  modulus,  and 
v*(p)  is  a  transform  domain  "Poisson's  ratio"  which  will  be  discussed 
later.  It  has  been  assumed  in  (436)  through  (438)  that  033  can  be 
neglected,  as  is  usually  done  in  theory  of  plates.  The  form  of  (436) 
through  (438)  follows  directly  from  the  discussion  given  in  Section  2, 
according  to  which  a  Laplace  transformed  viscoelastic  stress-strain  law 
is  completely  analogous  to  an  elastic  stress-strain  law  with  p- 
multiplied  transforms  of  relaxation  moduli  replacing  elastic  moduli; 
hence,  the  form  of  (436)  through  (438). 


♦ 


"u  =  p 


E  c  +  v  E  e 
a  11  a  1  22 


1  - 


T  -* 


*  Va 


22 


■  P 


=t  (c22  +  ”a  ,’n) 


l  -  -r 


*  Va 


°12  =  2p6a  ll2 


(439) 


(440) 


(441) 
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The  LT  of  the  bending  and  twisting  moments  in  the  plate  are 


given  by 


°ij  X3dx3  * 


(442) 


i.j  =  1.2 


Substituting  (435)  into  (439)  through  (441)  and  substituting  the  result¬ 
ing  expressions  into  (442),  we  find  that 


h3  ^  "'ll  *  \  ET  “'22 

12  s:  ,2 

1  -  ~rv 

E*  3 
a 


(443) 


M22  =  “  p  12  ~ 


h3  VVa  "»ll+"»22) 


i  Et*2 
a 


(444) 


h  '*  * 

M,  =  2  G  w,_„  . 

12  12  a  12 


(445) 


Here,  w,-^,  w>22  and  w»i2  denote  the  second  partial  derivatives  of  w 
with  respect  to  and  X2. 

The  quasi-static  equilibrium  equation  of  a  plate  in  terms  of 
bending  moments  is  given  by 


3\l  .  ..  s2h!2  .  *\l  . 

2  +  2  3x  3x  +  2  +  q(xl’  X2 *  *  °  * 

9x^  1  2  3x2 


(446) 


where  q  is  the  transverse  load  per  unit  area  of  plate,  positive  downward. 
See,  e.g.,  Timoshenko  and  Woinowsky-Kreiger  [16].  Taking  the  LT  of  (446) 
and  substituting  M  from  (443)  through  (445),  we  find  that 


**1111  4  +  2  **1122  2  2  +  **2222  4  =  q(xl’  x2*  p)  •  (447) 

9x^  9x2  9x2 


where 
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(448) 


(449) 


°2222 


(450) 


Equation  (449)  is  the  governing  differential  equation  of  the  LT  of  the 
deflection  of  the  plate. 

We  now  consider  a  typical  boundary  condition.  For  a  freely 
supported  rectangular  plate,  0  <  <  a,  0  <  <  b. 


w(x1>  0,  t)  =  w(xlt 


b,  t)  =  0 


(451) 


w,n(°,  x2,  t)  =•  (a ,  x2,  t)  =  w,22(x1,  0,  f)  =  w,22(xr  b,  t)  =  0  . 

(452) 

For  a  built-in  plate,  (452)  is  replaced  by 

w, 1 (0 ,  x2,  t)  =  w,1(a,  x2,  t)  =  w>2(x  ,  0,  t)  =  w,2(w1,  0,  t)  =  0  . 

(453) 

Since  our  governing  equation  is  for  the  LT  w  of  w,  the  bound¬ 
ary  conditions  (451)  through  (453)  have  to  be  Laplace  transformed. 

They  obviously  remain  in  the  same  form  after  Laplace  transformation. 

For  a  simply  supported  plate,  we  thus  have 

w(0,  x2,  p)  =  w(a,  x2,  p)  =  w(x  ,  0,  p)  =  w (x^ ,  b,  p)  «  0  ,  (454) 


w,j^(0,  x2 ,  p)  —  w,^(a,  x2 ,  p)  w»  (x^  0,  p)  —  w,22(x1,  b,  p)  =  0  . 

(455) 
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Suppose  now  that  the  load  on  the  plate  has  the  form 

q(x,  ,  x  ,  t)  =  q  (x,  ,  x9)  4>(t)  ,  (456) 

x  L  o-i .  i 

where  ( t )  is  any  function  of  time.  Then  the  Laplace  transform  of  (456) 
is  given  by 

q(x1,  x2,  p)  =  qQ  4>(p)  ,  (457) 

and  this  expression  now  has  to  be  used  on  the  right  side  of  (447). 

Let  qQ  in  (456)  be  expanded  into  the  double  Fourier  series 


qo<xl’  x2)  ■  l  £  aln  Vl  sl"  6nx2  1  (458) 

m^l  n=l 

where 


mfi 

m  a 


6  -  , 
n  b 


Tnn  ab 


q  (x. ,  x„)  sin  a  x.  sin  6  x„  dx,dx„  . 
^o  1  2  ml  n2  12 


a  solution  of  (447)  in  the  form 


00  00 


w(xr  x2 ’  P>  "  l  l  «mn  8in  Vl  8in  BnX2  ’ 

m=l  n*l 


(459) 

(460) 

(461) 


(462) 


It  is  seen  that  (462)  satisfies  all  the  boundary  conditions 

(454)  and  (455)  of  the  simply  supported  plate.  Now  substitute  (458) 

into  (457),  and  then  substitute  the  resulting  expression  and  (462)  into 

(447).  Equate  coefficient  of  sin  a  x.  sin  6  x„  on  both  sides.  This 
jju  ml  n2 

procedure  yields 


w 

mn 


(p) 


qmn  *(p) 


**  4  »*  2  2**  4 

Dllll  am  +  2D1122  am  6n  +  °2222  6n 


(463) 
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Equation  (463)  must  be  inverted  into  the  time  domain  to  yield 


w 

mn 


[w 

mn 


(p)] 


It  follows  from  (462)  and  (464)  that 


(464) 


w(x  ,  x- ,  t)  =  y  }  w  (t)  sin  a  x.  sin  B  x„  .  (465) 

1  z  lj1  L1  mn  ml  n  2 

m=l  n=l 

This  completes  the  formal  solution  of  the  problem  of  the  sim¬ 
ply  supported  rectangular  viscoelastic  plate. 

The  numerical  exploitation  of  the  solution  which  has  been 
given  here  unfortunately  involves  immense  difficulties.  The  quanti¬ 
ties  ^1122 »  an^  6*222  t*ie  denominator  °f  (463)  are  defined  in 

terms  of  transforms  of  viscoelastic  relaxation  moduli,  (448)  through 
(450)  ,  which  are  themselves  complicated  functions  of  transforms  of 
viscoelastic  relaxation  moduli  of  matrix  and  of  elastic  moduli  of  fi¬ 
bers.  Analytical  inversion  is,  in  general,  out  of  the  question,  and 
the  inversion  should  be  carried  our  numerically  by  the  use  of  compu¬ 
ters.  Here  no  attempt  will  be  made  to  carry  out  numerical  inversion  by 
computer.  Instead,  we  shall  approximate  the  viscoelastic  effective 
properties  of  the  fiber  reinforced  material  with  the  aim  of  making  an 
analytical  inversion  possible.  The  purpose  is  thus  not  to  solve  one 
case  precisely  from  a  numerical  point  of  view  but,  rather,  to  obtain  an 
estimate  of  the  viscoelastic  behavior  of  .  r  reinforced  plates  in 
bending.  We  shall  assume,  whenever  convenient,  that  the  matrix  is 
incompressible.  This  is  generally  a  fair  approximation  for  polymeric 
viscoelastic  materials. 

The  first  quantity  to  be  discussed  is  v*(p).  Its  interpreta¬ 
tion  is  as  follows:  suppose  that  the  fiber  reinforced  material  is  sub¬ 
jected  to  a  uniaxial  stress  constant  in  time  for  t>o,  a-jj(t)  =  o^  H(t) , 
in  fiber  direction.  This  produces  an  axial  average  strain  Eji(t)  with 
transform  fn(p)  antf  a  transverse  average  strain  e 2 2 ( t )  =  E33U)  with 
transform  £22(p)  =  ^33(p)*  transform  domain  Poisson's  ratio  va(p) 

is  defined  by 


-  * 


va(p) 


e22 ^P  ^ 

en(p) 


A  physical  time-dependent  Poisson's  ratio  may  be  defined  by 


(466) 
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(467) 


* 

v 

a 


(t)  =  - 


(t) 

To  • 


Note  that  v  (p)  is  not  the  transform  of  v  (t). 

It  has  been  shown  in  Reference  [3]  that  for  incompressible 
viscoelastic  matrix  and  stiff  elastic  fibers,  va(t)  is  time 
irdependent  and  is  given  by 


*  -  1  * 

v  (t)  =  —  v  +  v,v,  =  v  i'-  <  t  <  °°) .  (468) 

a  2  m  f  f  a  “ 

Here,  vm  and  v^  are  the  volume  fractions  of  matrix  and  fibers,  respec¬ 
tively,  and  vr  is  the  elastic  fibers'  Poisson's  ratio.  The  middle 
quantity  in  (468)  is  henceforth  denoted  v*. 

It  follows  from  (467)  and  (468)  that 

^22 =  ~  va  ?11^  ’  (469) 

Comparison  of  (469)  with  (466)  then  shows  that 

v  (p)  -  v  •  (470) 

a  a 


Next,  we  consider  Ea(p),  which  is  the  transform  of  the  axial 
Young's  relaxation  modulus.  It  has  been  shown  in  Reference  [3]  that 
E*  (t)  is  practically  time  independent  and  is  given  by 

E* (t)  =  <E  v  +  E  vf )  H(t)  =  E*  H(t)  ,  (471) 

a  m  m  r  t  a 

where  Em  and  E^  are  respectively  matrix  and  iber  elastic  Young's  moduli, 
and  Efl  without  time  argument  henceforth  denotes  the  middle  quantity  in 
(471).  Hence,  we  have  from  (471), 

l*(p)  -  f  .  (472) 

Now  the  transverse  Young's  modulus  Ej  of  a  uniaxially  fiber 
reinforced  material  is  much  smaller  than  the  axial  modulus  E*;  see 
Hashin  and  Rosen  [1].  The  same  is  true  to  a  larger  extent  for  a  visco¬ 
elastic  fiber  reinforced  material;  i.e., 
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>> 


(473) 


A 

E  (t) 
a 


E*(t) 


Equation  (473)  is  physically  obvious,  since  for  axial  loads  the  stiff 
elastic  fibers  take  the  major  part  of  the  load,  while  for  transverse 
loads  the  weak  viscoelastic  matrix  takes  the  major  part  of  the  load. 
Since 


0p) 


E  (t)  e  pt  dt 
a 


~  A 


et(p)  = 


E*(t)  e"pt  dt 


it  follows  from  (473)  through  (475)  that 


Efl(p)  >>  ET(p)  . 


A  typical  value  of  (468)  is  given  by 


(474) 


(475) 


(476) 


v  =  v  =  0.5,  v  =  0.2  . 
i  m  f 


Then 


v*  =  0.35  .  (477) 

a 

Using  (476)  and  (477)  in  (448)  through  (450),  it  is  seen  that  the  quan¬ 
tity  v*2  E^/E*  can  be  safely  neglected  in  comparison  to  unity.  Further¬ 
more,  we  introduce  (470)  and  (472)  into  (448)  through  (450).  We  then 
obtain  the  approximations 


=  TT  E 


1111  12 


°1122 


h  * 

U  >>  ‘“a  EI  +  2Ga>  ' 


2222 


12  P  ET 


(478) 

(479) 

(480) 
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We  now  have  to  consider  the  quantity  ££,  which  is  the  trans¬ 
form  of  the  transverse  Young's  relaxation  modulus  (normal  to  fibers). 
The  elastic  can  be  expressed  in  terms  of  other  moduli  in  the  follow¬ 
ing  way  (see  Reference  [2],  p.  16): 


4G„ 


*  ★  *  *2  * 

1  +  Gt/Kt  +  4G,rv//E 


T  a 


a 


(481) 


ft 

If  the  matrix  is  incompressible^or  even  nearly  incompressible,  then  Kj 
becomes  very  much  larger  then  G-j.,  and,  consequently ,  the  ratio  g£/k£  can 
be  neglected  in  the  denominator  of  (481).  In  addition,  for  stiff  fibers 
E*  is  generally  very  much  larger  than  G%  and  is  a  small  quantity. 

So  the  third  term  in  the  denominator  of  (481)  car  also  be  neglected.  We 
thus  obtain  the  approximation 


E*  -  4G*  ,  (482) 

valid  for  incompressible  matrix  and  stiff  fibers. 

In  transform  domain,  the  relation  (482)  becomes 

E*  -  4G*  .  (483) 

According  to  the  discussion  given  in  Section  8  for  incompressible  ma¬ 
trix  and  stiff  fibers,  the  transverse  relaxation  modulus  G*(t)  is 
approximately  given  by  (302) .  Taking  the  transform  of  this  equation  and 
introducing  it  into  (483),  we  obtain 


S  '  <  -  4  Ve)  •  <484> 

where  c  is  the  volume  fraction  of  fibers.  Consequently,  we  can  now 
further  approximate  (478)  through  (480),  obtaining 


z*  '  h3  * 
Dllll  12  Ea  ’ 


(485) 


D. 


122 


(1  +  2 v*)  , 

a 


(486) 


°2222 


(487) 
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We  now  assume  that  the  load  on  the  plate  is  constant  in  time 
for  t>o.  This  means  that  (J>(t)  in  (456)  is  given  by 


<Kt)  =  H(t)  . 


Consequently , 


4>  (p)  =  “  • 

P 

Introducing  (484)  through  (488)  into  (463)  ,  we  obtain 


(488) 


(489) 


12  Wmn(p) 


(-) 


*4  l  +  c~  *22  4  n 

E  a  +  4  G  (p)  [(1  +  2v  )  a  B  +  8  P 

am  l-cmr  amn  n 


(490) 


Equation  (490)  cannot  be  inverted  without  specification  of  the  func¬ 
tional  form  of  Gm(p).  We  make  the  simple  assumption  that  the  visco¬ 
elastic  matrix  is  characterized  by  a  simple  Maxwell  model.  In  that 
event,  p  Gm(p)  is  given  by  (286).  Introducing  this  result  into  (490), 
we  have,  after  rearrangement, 


where 


.3  q 

h  ,  .  mn 

12  Wmn(p)  =  TA 


1  +  T 


*  4  , 

mn  p  (E  a  /TA  +  p) 
v  a  mn  mn  r 

(no  sum  on  m,n) 


(491) 


*4  1  +  c  *22  4 

A  =  E  a  +  4G  7 - -  [1  +  2v  )  u  6  +  B  ] 

mn  am  ml-c  amn  n 


(492) 


Equation  (491)  is  easily  inverted  after  separation  of  the  right  side 
into  partial  fractions.  The  result  is 
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mn 
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mn 

h3/12 


r*  4 

i  E  a 

1  ,  a  m 
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mn 
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.  (493) 


Equations  (465)  and  (493)  specify  the  deflection  of  the  plate.  It  is  of 
interest  to  compute  the  initial  deflection  of  the  plate  w(x^,  x^,  o)  and 
the  final  deflection  w(x^,  >  °°)  •  From  (493),  we  have 
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mn 
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mn 
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h  E  ct  712 
a  m 


(495) 


It  is  seen  that  in  spite  of  the  fact  that  the  Maxwell  matrix  by  itself 
has  unbounded  creep  deformations,  the  plate  deformation  has  a 
finite  limit.  This  phenomenon  is  due  to  the  effect  of  the  fibers. 

Figure  24  shows  the  variation  with  nondimensional  time  t/T 
of  the  deflection  at  the  center  of  a  square,  simply  supported,  visco¬ 
elastic  fiber  reinforced  plate.  The  plate  and  material  data  are  as 
folJows : 

E  =  5.0  x  10^  psi 
m  v 


Ef  =  10.5  x  106  psi 
vf  =■  0.2 


G  =  1.7  x  105  psi 
m 

c  =  0.5 
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Figure  24.  Late 


13.  AX  I  SYMMETRICALLY  LOADED  VISCOELASTIC 
FIBER  REINFORCED  CIRCULAR  CYLINDRICAL  ShELL 


The  shell  is  shown  in  Figure  25  with  its  coordinate  system. 
Axisymmetry  under  axisymmetric  loading  requires  that  material  properties 
also  be  axisymmetric.  For  uniaxial  reinforcement,  this  is  possible  only 
for  two  cases:  when  the  fibers  are  in  generator  direction  or  when  the 
fibers  are  in  circumferential  direction,  normal  to  generators.  In  view 
of  previous  notation  for  material  properties  where  x^  always  was  in 
fiber  direction,  we  can  here  make  the  following  identification: 

Longitudinal  reinforcement  x  -*■  x^  0  ■+  X2  (496 ) 

Circumferential  reinforcement  x  -*■  0  -*  x^  .  (497 ) 


We  consider  a  static  shell  whose  only  load  per  unit  shell  area 
is  in  z-direction,  normal  to  the  shell  surface, 


q2  “  q  =  q(x)  . 


(498) 


In  the  axisymmetric  case  the  shear  membrane  force  Nx0 ,  the  twisting  mo¬ 
ment  Mxe,  and  the  shear  force  Q0  all  vanish.  Furthermore  there  is  no  0 
dependence  in  any  quantity.  Simplifying  the  shell  equilibrium  equations 
accordingly  (see.e.g.,  Reference  [17]),  we  obtain 


dN 
_ x 

dx 


o 


(499) 


dQ  N 

dT  "  T  +  q(x)  "  0  *  (500> 


dM 

dT  '  Qx  "  0  ’  <501> 

where  Nx,  Ne  are  normal  membrane  forces,  Mx  is  normal  bending  moment, 
and  Qx  is  the  shear  force.  Elimination  of  Qx  from  (500)  and  (501)  yields 

A  N0 

— r  “  ir  +  q(x)  "  0  »  (302> 

dx 

while  (499)  implies  that 


N  =  N  =  const, 
x  0 


(503) 
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Figure  25.  Cylindrical  Shell  Coordinate  System. 


In  order  to  simplify  the  discussion,  we  shall  derive  the 
equations  of  an  elastic  fiber  reinforced  shell,  and  we  shall  then  proceed 
to  the  viscoelastic  case  via  the  correspondence  principle.  On  the  basis 
of  the  Kirchhoff  assumption,  axisymmetry,  and  the  thin  shell  approxi¬ 
mation  h/R  <<  1,  the  only  surviving  strains  in  the  shell  are 


du  d^w 

Cxx  dx  “  z  ,  2  » 
dx 


(504) 


E 


ee 


(505) 


Here,  z  is  measured  from  the  shell  middle  surface,  normal  to 
it  and  positive  outward,  u  and  w  are  displacements  in  x-  and 
z-directions,  respectively. 


Let  it  now  be  assumed  that  the  reinforcement  is  longitudinal. 
Then  (496)  applies,  and  the  stress-strain  law  is  given  by 
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*  * 

where  Ea  is  the  axial  Young's  modulus,  va  is  the  axial  Poisson's  ratio, 
and  Ej  is  the  transverse  Young's  modulus.  Solving  (506)  and  (507)  for 
the  stresses  yields 
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where 
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(512) 


The  shell  membrane  forces  and  bending  moments  are  given  by 
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(513) 

(514) 


We  now  introduce  (504)  and  (505)  into  (508)  and  (509),  and  we  introduce 
the  resulting  expressions  into  (513)  and  (514).  Thus,  shell  membrane 
forces  and  bending  moments  are  obtained  in  terms  of  u  and  w.  These  are 
in  turn  introduced  into  the  equilibrium  Equations  (502)  and  (503).  This 
procedure  leads  to 
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Elimination  of  yields  the  equation 
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(518) 


The  membrane  force  NQ  in  (517)  is  equal  to  the  constant  pre¬ 
scribed  membrane  force  In  x-directions  at  the  end  sections  x  *  0,  L  of 
the  cylinder.  If  there  is  no  membrane  force  prescribed,  NQ  vanishes. 

Usual  boundary  conditions  at  x  •  0,  L  prescribe  any  two  of  the 
dw 

quantities  w,  Mx,  and  Qx< 

*  * 

For  a  viscoelastic  fiber  reinforced  shell  cylinder,  Ea,  ET, 
and  v*  have  to  be  replaced  by  p^a(p),  pEj(p),  a°d  va(p)  an  analogy  to 
the  treatment  given  in  Section  12.  The  load  q  is  now  q(x,t),  a  func¬ 
tion  of  time  also,  and  has  to  be  replaced  by  its  transform  q(x,p). 

Let  it  be  assumed  for  simplicity  that  the  cylinder  is  not  sub¬ 
jected  to  longitudinal  membrane  load  at  the  ends.  Then  (517),  written 
for  the  Laplace  transform  of  w(x,t)  which  is  denoted  by  w,  becomes 


.4* 

d  w  .  - ,  x* 
— 7  +  y(p)w 

dx 


i2fL 


h3  A  (p) 
xx  r 


(519) 


where 


y(p) 


12  \«(P>  *88  W  ~ 

R2h2  A*  <p) 


(520) 
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i 


and 


PE„ 


A  (p) 
xx  r 


AX0(P) 


A*  * 

E  .2 

T  -* 

1-?v* 

a 


pv«  et 

A* 

,  fi-*3 

E*  V* 

a 


(521) 


(522) 


*ee(p) 


(523) 


1  -  -r  v 


For  very  stiff  fibers  and  nearly  incompressible  viscoelastic 
matrix,  the  approximations  used  in  Section  12  apply  in  a  completely 
analogous  way.  These  are  contained  in  Equations  (4*'e>  through  (487)  for 
the  approximate  TD  plate  stiffnesses.  In  the  present  case,  Equations 
(521)  through  (523)  become 


*  ,  *  * 

A  (P)  -  E  » 


xx 


*  *  *  A* 

A  Q(p)  ■  4v  p  G  , 
x9  r  ar  a  * 


A00^ 


4pG„ 


(524) 


(525) 


(526) 


Introduction  of  (524)  and  (526)  into  (520)  yields 

48p  8* 

Y(P)  -  Tji  . 

R  ti  E 

a 


(527) 


In  the  elastic  case,  Equation  (517)  is  easily  solved  in  gen¬ 
eral,  and  the  solution  is  then  fitted  to  the  boundary  conditions.  It 
should,  however,  be  borne  in  mi. id  that  the  homogeneous  solution  of  (517) 


.f 
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1/4 

contains  ihe^terms  exp[(-y)  ].  In  the  viscoelastic  case,  y  has  to  be 
replaced  by  y(p)  as  given  by  (527),  Thus  the  Laplace  transform  w  will 
be  obtained  in  a  very  complicated  form,  tnd  the  inversion  will  prove  to 
be  extremely  difficult.  Indeed,  analytical  inversion  is  out  of  the 
question.  It  is  therefore  more  convenient  to  use  a  Fourier  method  for 
the  solution  of  (519),  since  this  results  in  a  much  simpler  form  of  w. 

Let 


q(x,t)  -  l  q  (t)  sin  ax  , 
-  n  n 

n-l 


n* 


n  L  * 


Then 


q(x.p)  ■  I  q„(p)  Sin  ax 
,  n  n 

n-l 

Assume  that  a  solution  of  (519)  is  the  form 


a  r 

w  -  )  w  (p)  sin  a  x  . 

,  n  n 

n-l 


(528) 

(529) 


(530) 


(531) 


Introduce  (530)  and  (531)  into  (519),  also  taking  into  account  (524). 
Equating  terms  of  identical  sine  functions,  we  have 


12  qr 


wn(p)  -  7T7 


h^E  Iy(p)  +  <\*J 

41  It 


(532) 


where  y(p)  is  given  by  (527).  The  solution  in  form  (531)  satisfies 
boundary  conditions  of  simple  support  at  x  -  0,  L;  i.e., 


w(x,t)  -  o 


d 

dx2 


A 

w(x,p)  -  o 

d.2y.foiP)  . 

dx2 


)  X 


0,  L. 


(533) 


Let 
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(534) 


w  (t)  -  LT"1  w  (p). 
n  n  r 

Then  the  solution  to  the  problem  of  the  simply  supported  shell  is  given 
by 

oo 

w(x,t)  ■  l  w  (t)  sin  a  x  .  (535) 

-  n  n 

n-1 

Suppose  now  that  the  loading  is  uniform  and  does  not  vary  in  time.  Then 

q(x,t)  -  qQ  H(t)  , 

q(x,p)  *  q0  “  •  (536) 

The  Fourier  expansion  of  1  in  o  i  x  <  L  is  given  by 


1  ■  —  5"  “*  sin  a  x  . 

x  S  i  -  n  n 
n-1, 3, 5 


(537) 


Hence,  in  the  present  case, 

4q 

_jo 

qn  imp 


n  -  1,  3,  5  —  . 


(538) 


Suppose  again  that  the  matrix  is  represented  by  a  simple  Maxwell  model. 
Then  as  in  Section  12, 


r*  1  +  c  nm^ 

pGa  "  1  -  c  1  +  Tp  * 


(539) 


Introducing  (539)  into  (527)  and  introducing  the  resulting  expression  to¬ 
gether  with  (538)  into  (532),  we  obtain 


wn(p) 


12 _ (1  +  Ip)2 _ 

*  3  2  2  4  2  * 

E  h  p{$Tp  +  $BT  p  +  a  (1  +  Tp)  } 
s  n 


(540) 


where 


4> 


12  4  (1  +  c) 

R2h2  1'c 


(541) 
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(542) 


8 


1  +  c  m 


a 


Here,  Gm  is  the  initial  matrix  shear  modulus.  Equation  (540)  can  be 
written  in  the  form 


o  (p)  ■  -U_  .1.  +  2Ip,  *  *V 

nvp'  *  3  4  2 

Eah  p  )i  +  isa  +  £_  ) 


u 


(543) 


where 


—  “Til 


+  (a) 


(544) 


2  +  4- 

Z  4 
a 

n 


2/  1+4 

4 
a 

n 


The  Laplace  inverse  transform  of  (543)  is 


(b) 


...  12  “>o 

wn(t)  '  TO  4 

E  h  Tina 
a  n 


1  + 


A7 


((1  +  2T;o>1  -  T2ui2  +  2T2c2u)2)2 


+  u)2  (1  -  c2)  (2Tcw)2]172  e_i>t  sin  (u^  Jl  -  *  t  +  *) 


(545) 


where 


if/  “  tan 


-1  2Toj</  1  -  t2 


5.  -U  - T^) 


T2oo2  (2c2  -  1)  +  1  -  2TCu 


-  tan 


(5W 


-  c 
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